PB87-181483 


Testing and Validating the CERES-Wheat 
(Crop Estimation through Resource and 
Environment Synthesis-Wheat) Model in 
Diverse Environments 


(U.S.) National Aeronautics and Space 
Administration, Houston, TX 


Dec 86 






INBHH 





REPORT DOCUMENTATION *• *> *• ' 

PAGE YM-15-00407-JSC-20244 

3. Aeciplanfa A creation No. 

mr 18148 8MS 

4. TWa and Subtitle 

Testing and Validating the CERES-Wheat Model in Diverse 
Environments 

S. Report Dote 

December 1986 

s. 

7. Authorial 

S. Otter-Nacke, D. C. Godwin and J. T. Ritchie 

a * - « - - ^ » — ^t — u— 

fwiipC. pvo* 

YM- 1 5-00407- JSC-20244 

L Mormlni Organization Noma and Addrasa 

US Department of Agriculture-Agriculture Research Service 
Grassland, Soil and Water Research Laboratory 
P.O. Box 6112, Temple, TX 76503-6112 

10. Proiect/Taek/WMi Unit No. 

NA 

11. ContmctfO or Grant(G) No. 

(O 

_ NA 

«5) 

12- Spon wring Orgamufon Name and Addraas 

AgRISTARS 

USDA-ARS 

13. Typo of Report O Period Covered 

Technical Report 
1980-1985 

14. 


IS. Su ppremsn ts ry Notes 


IS. Abstract (Limit: 200 srorda) 

CERES-Wheat is a computer simulation model of the growth, development, and 
yield of spring and winter wheat. It was designed to be used in any location 
throughout the world where wheat can be grown. The model is written in FORTRAN 77, 
operates on a daily time stop, and runs on a range of computer systems from micro- 
computers to mainframes. Two versions of the model were developed: one, CERES- 

Wheat, assumes nitrogen to be nonlimiting; in the other, CERES-Wheat-N, the effects 
of nitrogen deficiency are simulated. This report provides the comparisons -f 
simulations and measurements of about 350 wheat data sets collected from throughout 
the world. Comparisons are made of model performance to predict grain yields, 
yield components, leaf area index and biomass, phenology, soil water balance, soil 
nitrogen balance and nitrogen uptake by plants. Agreement between model output 
and experimental results were acceptable in most instances. Statistical analysis 
indicated that some problems have yet to be resolved with the prediction of N 
uptake and grain protein concentration in some data sets. 


17. Document Analysis a. D es c ript o rs 

computer simulation model, wheat, nitrogen deficiency, grain yields, 
yield components, leaf area index, biomass, phenology, soil water balance, 
water deficiency 

b- lds n t H lai t /Opeo<ndtd Tams 


e. COSATI FMd/Oroup 


1SL Availability Statement 

19. Security Claes (This Report) 

21. No. of Pagoa 

Release unlimited 

Unclassified 

155- 


20. Security Claaa (This Page) 

22. Price 


Unclassified 



(See ANSI-Z39.lt) Sat Inst ructions on Ravers* OPTIONAL FOMI 272 (4-77) 

(Formerly NTIS-3S) 
Papartmant of Commar c a 


TESTING AND VALIDATING THE CERES-WHEAT MODEL 


IN DIVERSE ENVIRONMENTS 


Testing and Validating the CERES-Wheat Model 
in Diverse Environments 


S. Otter-Nacke-*- 
D . C . Godwin^ 
J. T. Ritchie 3 


DSDA-ARS 

Grassland, Soil and Water Research Laboratory 
P.0. Box 6112 
Temple, Texas 76503-6112 
USA 


International Fertilizer Development Center 
P.0. Box 2040 

Muscle Shoals, Alabama 35662-2040 
USA 


Post-Doctoral Fallow, Department of Crop and Soil Sciences, 286 Plant and Soil 
Sciences Buildlnt, Michigan State University, East Lansing, MI 18824; formerly Research 
Associate, Texas Agricultural Experiment Station, P.O. Box 6112, Tomple, TX 76503. 

2 

Agronomist/Systems Modeler, Agro-Economic Division, International Fertilizer 
Development Center, P.O. Box Z040, Muscle Shoals, AL 35662. 

^Professor, Homer Nowlin Chair, Department of Crop and Soil Sciences, 286 Plant and 
Soil Sciences Building, Michigan State University, East Lansing, MI 48824; formerly 
Research Leader, Crop Systems Evaluation Research Unit, USDA-ARS, P.O. Box 6112, Temple, TX 
76503 . 


• i • 

III 



CONTENTS 


INTRODUCTION 1 

FEATURES OF CERES -WHEAT 1 

CERES -WHEAT (NON-NITROGEN VERSION) 2 

Data Base 3 

Methods 7 

Results and Discussion 8 

Summary 27 

CERES -WHEAT NITROGEN 27 

Data Base 28 

Individual Data Sets 29 

Results and Discussion 32 

Summary 46 

APPENDIX A 47 

CERES -Wheat Model Validation Results From Individual 
Experiments 48 

APPENDIX B 106 

CERES -Wheat-N Model Validation Results From Individual 
Experiments 107 

BIBLIOGRAPHY 144 


v 


Precede page blank 



ILLUSTRATIONS 


Figures 

1. Location of Experimental Sites 4 

2. Predicted Versus Measured Yield 16 

3. Predicted Versus Measured Grain Weight 16 

4. Predicted Versus Measured Number of Grains 17 

5 . Predicted Versus Measured Number of Ears 17 

6. Predicted Versus Measured Above Ground Dry Matter 18 

7. Predicted Versus Measured Leaf Area Index 18 

8. Comparison of Predictions of the CERES -Wheat-N 

Model with Observed Data from Field Experiments 38 

Tables 


1 . 

Data Source for Model Testing 


5 

2. 

Summary Measures for Independent Data 

Sets 

19 

3. 

Summary Measures for Dependent Data Sets 

19 

4. 

Difference Measures for Independent Data Sets 

19 

5. 

Difference Measures for Dependent Data 

Sets 

19 

6. 

Listing of Final Results of Independent Data Sets 

20 

7. 

Listing of Final Results of Dependent 

Data Sets 

24 

8. 

Summary Measures for all Data Sets 


40 

9. 

Effect of Errors in N Uptake on Errors 
Apparent Recovery 

and 

35 

10. 

Listing of N-Model Testing Data Base 


41 


vii 


Preceding page blank 



ACKNOWLEDGEMENTS 


The authors wish to thank the numerous scientists, their 
collaborators, and institutions whose data enabled the 
development and validation of the models. Primary data was 
obtained from the following institutions: 


Commonwealth Scientific and Industrial Research Organization 
(CSIRO) - Australia 

New South Wales Department of Agriculture - Australia 
Western Australia Department of Agriculture - Australia 
Agriculture Canada - Canada 

Agriculture and Food Research Council - England 
Rothamsted Experiment Station - England 
University of Nottingham - England 

Institut National des Recherches Agronomiques (INRA) - France 
Kibbutz Sole Boker, Negev Desert - Israel 

International Maize and Wheat Improvement Center (CIMMYT) - 
Mexico 

Central Agricultural and Biological Organization (CABO)- 
Netherlands 

Soil and Irrigation Research Institute - South Africa 
International Center for Agricultural Research and Dry Areas 
(ICARDA) - Syria 

Agricultural Research Services -United States Department of 
Agriculture (ARS/USDA) - USA 
Kansas State University - USA 
Montana State University - USA 
North Dakota State University - USA 
Oregon State University - USA 
Texas A&M University - USA 
Washington State University - USA 
Technische Universitat, Mvinchen - West Germany 
Deutscher Wetterdienst Agrarmeteorologische Beratungs und 
Forschungsstelle Ahrensburg - West Germany 


We would also like to acknowledge the assistance of Reah 
Janise Battenfield in editing this manuscript, Terese 
Heineman Baker for its preparation, and Mat Krogulecki and 
Tom McCafferty for preparing the computer -generated graphs. 


lx 


Preceding page blank 


ERRATA PB181483; Testing and Validating the CERES -Wheat 
Model Diverse Environments. 


To be added after last paragraph of page 31 and before the 
first paragraph of page 32. 


Dufur. Pendleton, and Madras Oregon. 1971 

These experiments were reported by Ambler (1974) and 
involved the comparison of five varieties of winter wheat 
over several different rates of fertilizer at several 

different sites in Oregon. Data from sites where climatic 
data were incomplete were rejected. The three remaining 
sites differ in annual rainfall and altitude. Solar 

radiation data from Klamath Falls were used for the Madras 
data set. Few significant variety x N interactions were 
recorded and so data from only two cultivars (Hyslop and 
Nugaines) were utilized. At Madras four rates of N (0, 90, 
180, and 270 kg N/ha) at planting with a second application 
of 90 kg N/ha applied at the booting stage was investigated. 
Initial levels of nitrate in the profile were high and the 

response to N was small. The experiment was irrigated. At 

the Dufur site, three rates of N (0, 17, and 34 kg N/ha) were 
used and at the Pendleton site four rates of N (0, 34, 67, 
and 101 kg N/ha) were used. No significant response to N was 
recorded at either site. 

Waite Inscitute. 1958 

These data were reported by Barley and Naidu (1962). The 
experiment examined the response to N of two varieties (Gabo 
and Bencubbin) . The rates of N used were (0, 33, 67, and 174 
kg N/ha) as ammonium sulfate. Solar radiation data were 
estimated from the recorded values of hours of sunshine. The 
soil was reported as a red-brown earth (Alfisol) . 

Bozeman. Montana. 1977 

These data were reported by Christianson and Killorn 
(1981) . This experiment examined differences in fertilizer 
use efficiency for applications made at different times. 
Nitrogen was applied as ammonium nitrate after seeding or 
broadcast several hours prior to a sprinkler irrigation to 
simulate application of fertilizer through the sprinkler 
system. The study investigated the effects of four rates of 
N applied at planting (0, 50, 100, and 150 kg N/ha) and five 
split-application patterns with a total application rate of 
either 100 or 125 kg N/ha. The soil was a deep silt loam and 
N responses were very marked. 
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For the independent data sets, an "index of agreement" 
(d) was calculated according to the procedure of Willmott 
(1982). This calculation was only performed on yield data. 
All statistical computations were performed with the 
statistical package SAS (1982). 

The second and much more extensive part of the results 
are graphics, showing the proximity of the time course of 
model output and experimental measurements. Unfortunately, 
most data sets provided no information available on the 
variance of published measurements. This would be of certain 
interest because observed values are not true values in a 
statistical sense, but are mean values of a very small sample 
with more or less large intrinsic error. 


Results and Discussion 

Summary Statistics . Tables 2 and 3 are summa zed 
measures for the dates of anthesis and physiolog. cal 
maturity, yield, grain weight, grain numbers per m^ (GPSM) , 
ear numbers per m , above ground dry matter, and leaf area 
index (LAI). Observed and predicted means are listed to show 
the difference between them. Standard deviations of observed 
and predicted values are listed to be checked for similar 
ranges. The intercept and slope of the regression equation 
with 0 as an independent variable provides information about 
the accuracy of predictions. Tables 4 and 5 contain the 
difference measures for the same list of variables. 

Difference Measures . While summary measures try to 
describe the quality of simulation, difference measures cry 
to locate and quantify errors. Difference measures fot the 
same variables as above are listed in Tables 4 and 5. The 
mean bias error (MBE) is regularly considerably smaller than 
the mean absolute error. The two values, being close, 
indicate the prediction to be biased in one direction, as is 
true for ear number and dry matter predictions with the 
independent sets . 

A negative MBE occurs when predictions are smaller in 
value than observations. Considerable errors of this kind 
are found with GPSM, ear number, and dry matter predictions 
of the dependent data sets, and with grain weight and date of 
anthesis of the independent sets. For all others, 
predictions are larger than observed values. 
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Since the value of the liBE is related to the range of 
values, it is not reasonable to compare across variables. 
However, comparisons between the two groups of data sets 
might show the necessity of requiring inde' -ndent data sets 
if it could be shown that the MBE and MAE ot independent data 
were always larger. Unfortunately, while true for yield, it 
does not apply for all variables; maximum error in predicting 
yields is smaller in the independent set. 

Willmott (1982) indicated that the closer the root mean 
square error (RMSE) to 0, the better the model. Certainly 
this applies to the systematic error. The unsystematic RMSE 
should approach RMSE in the system: 

MSE - MSE s + MSEu 

The tendency of RMSE s toward 0 is obvious, compared to RMSE. 
Since the mean is involved in the computation of this 
measure, comparisons across groups of data sets are not 
advisable. 

The values of RMSE s and RMSE, computed to predict yields 
for the dependent data set, aie very similar: 25377 and 
25424, respectively. Thus, compared to RMSEu of 1552, the 
model meets the above requirement. 

Willmott (1982) also suggests the computation of d, a 
quality me ure for models, or index of agreement: 

d - 1 - [2 (F.- 0.) 2 / 2 (|P'| + | 0 ' | ) 2 ] , 0 < d < 1 
i-1 1 i-1 

where - P^ - 0 and 0^ - 0^ - 0 

For a good model, d should approach 1. For CERES -Wheat yield 
predictions, d equals 0.8825. This quality measure is mainly 
used to compare different models, but it provides additional 
useful information of model performance when used to evaluate 
a single model. 

Descriptive Statistics . Descriptive statistics based on 
the values shown in Tables 6 and 7 are reported separately 
for the two groups of data sets. 




ERRATA PB181483 ; Testing and Validating the CERES -wheat Model in Diverse 
Environments . 

This pace replaces page 10 of the original document . 


Phasic Development . The accuracy of simulating the 
phasic development of a crop is crucial fcr the accurate 
simulation of crop growth and yield. Therefore, this part of 
the simulation received considerable attention. Figures 2-7 
show the scatter of predicted versus measured variables, 
including the 1:1 line +1 standard deviation and the matching 
regression. 

Regression equations for anthesis and physiological 
maturity show an excellent fit with slopes close to 1 and 
high regression coefficients. Observed means are 
considerably different for the two groups of data sets, the 
independent data sets reading anthesis and maturity 12 and 19 
days, respectively, later. This reflects the fact that the 
origin of this subset of the data base is relatively far 
north, where the wheat crop season is extended into the 
latter part of the year. 

Difference measures ir. Tables 4 and 5 indicate small 
mean bias errors and mean absolute errors of five and seven 
days, respectively. The maximum errors for the two dates in 
the dependent subset occurred in the same run; therefore, the 
time available for grain filling was correct. When the date 
of anthesis was predicted too early, maturity came too early 
also. 


Yield. Grain Weight, and Grain Numbers/m ^ . The 
independent data set is characterized by a higher mean yield 
caused by higher grain weight and slightly bigger grain 
numbers. This confirms that the data sets came from 
different locations. High yields were mostly from northern 
wheat growing areas. However, observed and predicted mean 
yields are very close resulting in small mean bias errors. 
The mean absolute error is in the range of 22 percent of the 
observed mean yield and GPSM, but only 15 percent of the 
observed grain weight. 

Figures 2-4 demonstrate the scatter of predicted versus 
measured data points around the regression line, the 1:1 line 
and within the limits of +1 standard deviation of the 
observed mean. The percentage of data points outside of 
these limits is small as to the total number of observations. 
This demonstrates the model’s reliability in simulating 
yield. Comparison of the regression to the 1:1 line proves 
the good fit, especially for the independent data sets. The 
regression coefficient, however, does not indicate a 
significant superiority of yield predictions over the rest of 
the data base (0.633 compared to 0.617). 
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It is interesting to note that measured and predicted 
mean yields within a subset did not differ nearly as much as 
between the two subsets. 

Final Ear Number. Final Dry Matter, and Maximum LAI . 
Additional variables to be checked are final ear number, 
final dry matter, and maximum LAI observed during the growing 
season. Since the average experiment had about five 
measurements in the season, there is no guarantee that a 
measurement was taken at the peak of leaf area development. 
There are relatively few observations in total. Therefore, 
observed maximum LAI has to be considered with precaution. 

Despite these restrictions, the mean bias error is not 
very large. The mean absolute ^rror&, however, amount to 40 
percent and 95 percent of the observed mean. Also, the 
regression coefficient is not very high (Figure 7) . (Please 
refer to the section where predicted and observed LAI are 
compared over the course of time for particular experiments 
to demonstrate the model's capability to simulate the canopy 
development . ) 

Leaves comprise most of the total amount of above ground 
dry matter (Figure 6) . Since final dry matter is highly 
correlated to yield, it must be accurately simulated. The 
model overpredicts dry matter (Table 4) by about 380 g/m^, 
which is considerably less than 1 standard deviation (525.5 
g/m^). In the regression, the model accounts for 41.1 
percent of the variation in the real world. 

For the group of independent data, simulation of final 
ear number (Figure 5) resulted in a very small slope (0.1398) 
and corresponding coefficient of determination (r^ - 0.014). 
The regression coefficient for the dependent data sets is 
0.40 and thus is higher than for dry matter. 

Measured and Model Estimates of Cron Growth . To ensure 
accuracy, any simulation of plant parts should be checked 
intensively against real world observations. The appendix 
contains comparisons between model output with all available 
experimental measurements. For some, real world data .ere 
obtained by digitizing graphs in publications, which adds 
uncertainty to the reliability of these data. These 
comparisons should be viewed with a critical mind. Figures 
in Appendix A are in the order of locations. 
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One problem with the data is that often only a few 
samples were taken during the growing season and for typical 
optimum curves like LAI or tiller number, the peak may be 
missing (see tiller number for Rutherglen, Australia) . The 
experiment in Rutherglen is one of the few wher n root 
measurements were taken. It demonstrates good correspondence 
of compared data despite high variations in measurements and 
a slightly too high growth rate in the model run. 
Partitioning was correct for this particular experiment as 
shown through a large number of samplings of root, stem, and 
ear weight. One discrepancy is that stem weight measurements 
seem to include leaves or sheaths, whereas the graph of the 
simulation shows stem weight only. 

Five, well-timed samples during the growing season 
demonstrated agreement with the tillering pattern at 
Murumbateman, Australia, in 1977. Biomass and LAI were 
simulated with high accuracy. The 1980 run at the same 
location did not perform quite as well, although the same 



TESTING AND VALIDATING THE 
CERES -WHEAT MODEL IN DIVERSE ENVIRONMENTS 

by 

S. Otter-Nacke, D. C. Godwin and J. T. Ritchie 


INTRODUCTION 

CERES-Wheat^- is a computer simulation model of the 
growth, development, and yield of spring and winter wheat. 
It was developed through the assistance of several scientists 
from the United States and other countries. CERES, Crop 
Estimation through Resource and Environment Synthesis, is 
derived from the Latin word for cereal. 

The model is written in FORTRAN 77, operates on a daily 
time step, and runs on a range of computer systems from 
micro -computers to mainframes (Ritchie and Otter 1985) . Two 
versions of the model were developed: one, CERES-Wheat, 
assumes nitrogen to be nonlimiting; in the other, CERES- 
Wheat-N, the effects of nitrogen deficiency are simulated. 
Although the models had been intensively tested and 
calibrated, they had not been validated. Before any model 
can be used with confidence, adequate validation or assess- 
ment of the magnitude of the errors that may result from its 
use should be performed. Model validation in its simplest 
form involves a comparison between simulated values and real 
world values. This paper describes the validation of both 
models and is presented in two sections. The first section 
(pages 2-46) summarizes the data from all the experiments. 
The last section (Appendix A and B) provides detailed 
comparisons from individual experiments . 


FEATURES OF CERES-WHEAT 

CERES-Wheat was designed to be used in any location 
throughout the world where wheat can be grown. It describes 
the major environmental and edaphic factors affecting yield 
by way of s Emulating the following processes : 

Phasic development or duration of growth stages on 
a thermal time scale as related to plant genetics, 
weather, and other environmental factors. 


^Documentation of the CERES-Wheat model is currently in 
review and will be published in the same series as CERES - 
Maize, A Simulation Model of Maize Growth and Development , 

C. A. Jones and J. R. Kiniry, Eds. 
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- Apical development as related to morphogenesis of 
vegetative and reproductive structures. 

Extension growth of leaves and stems and 
senescence of leaves . 

Assimilate accumulation and partitioning. 

- Soil water status and its effect on growth and 
development . 

Nitrogen status and its impact on growth and 
development . 


GERES -WHEAT (NON-NITROGEN VERSION) 

Calibration, validation, and sensitivity analysis are 
three different ways of evaluating models. Logically, 
validation is the last step in model testing. Although 
important to proper model evaluation, calibration cannot be 
the only or last step. Penning de Vries (de Vries and van 
Laar 1982) stated that calibration is a "very restricted form 
of evaluation," an "adjustment of some parameters such that 
model behavior matches one set of real world data," and that 
extensive calibration "degrades simulation into curve 
fitting." 

Sensit* vity analysis reflects the effects of perturbed 
model pa* meters on model output, such as biomass and yield. 
The information can be used to determine whether the model is 
overly sensitive or not sensitive enough to certain processes 
being simulated. 

Final validation of a model can only be conducted after 
a. calibration and a sensitivity analysis. Validation 
requires independent data sets to verify the behavior of the 
model under real world conditions. In its simplest form, it 
is ,1 comparison between simulation and the real world. 

statistical methods can provide some useful information 
on model performance. The correlation coefficient (r) and 
coefficient of determination (r^) are of limited value. 
•Jillmott (1982) contends that r and r^ are of little 
practical value in evaluating model performance, because 
their magnitudes are not consistently related to the accuracy 
of prediction. Better criteria are bias, mean bias error, 
varianc*. of the distribution of differences, root mean square 
erro v (RMSE) , and mean absolute error (MAE) . He recommends , 
ar minimum, reporting predicted mean (P) and observed mean 
(0) , the standard deviation of the predicted (Sp) and 
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observed (S Q ) variable, the intercept (a) and slope (b) of 
the least-squares regression - a + b 0^ (Willmott 1982). 

Eight criteria were established by AgRISTARS for testing 
and evaluating crop yield models: mean square error, 
variance bias, proportion of years beyond a critical error 
limit, worst and second to worst performance during the 
testing period, range of accuracy, direction of change from 
mean yields, and correlation coefficient between actual and 
model predicted yields for a set of independent test y>ars 
(Wilson and Sebaugh 1981) . Other criteria include mean 
absolute difference and length of an 80% prediction interval 
for an observation of the difference between observed and 
calculated yield (Kornher and Torsell 1983). When defining 
the limits of acceptance for model performance, McMahon 
(1983) claims that 85S of the variance in crop yield should 
be accounted for or the predicted mean should be within one 
standard deviation of the observed mean. 

A temptation when calibrating the model is to use some 
of the information from data sets collected for validation. 
Results of these tests are of limited value, because experi- 
mental data may have been used to evaluate one particular 
equation in the source code. For CERES -Wheat, other data 
sets were run during the model's development. Thus, model 
weaknesses were based on more than one experiment, which 
allowed improvements to be made, and the same data base could 
be used for further validations. In the validation process 
discussed in this paper, the data base is divided into (1) 
independent data sets and (2) previously used data bases. 


Data Base 

A good data set should contain as much detail as 
possible to describe the process of plant growth. However, 
taking a sufficient number of measurements throughout the 
season can make a simple experiment expensive and location- 
specific. A universal model must be tested at many 
locations . 

To test the CERES -Wheat model in the northern and 
southern hemisphere, a data base was assembled to represent a 
diversity of environments, including short growing season 
spring wheat crops, environments with limited water avail- 
ability, subtropical wheat growing areas with little vernali- 
zation, and regions with temperature extremes. The data base 
of 283 data sets (collected between 1960 and 1983, from 25 
different sites , includes published experimental data from 
the United States as well as from other countries as far 
south as Australia and as far north as West Germany, repre- 
senting a range of latitudes from 36'S to 54°N (Figure 1). 
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Figure 1. Location of Experimental Sites. 




Validation depends on the quality of data available. To 
provide a high-quality data set for wheat model testing, 
extensive crop measurements throughout the season are 
required. As most research work is for purposes other than 
collecting data for testing a model, published data are not 
as comprehensive as desired. Thus, in many instances, some 
compromises for model input were often necessary to fully 
establish the data base. Minimum input included daily 
weather data for the entire season, soil water information, 
phenological observations, yield, and yield components 
(Table 1). 

Since several data sets in the data base had been used 
to calibrate and validate former versions of the CERES model, 
and statisticians emphasize that only truly independent data 
should be used for validation, the data base was divided into 
two groups: 

1. Dependent data sets (i.e. data sets used during 
model development) (130) . Some data were used 
to calibrate parts of a particular subroutine of 
the model. Results of their testing will be 
shown, however, because other parts of the model 
may be unaffected despite some modifications 
from that experiment. 
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Table 1. Data Sourea for Modal Tearing. 
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Avignon ' 79 n 

2 

2 




** 

** 


Bolgnevllle '76 n 

4 

2 



2 

*** 

** 


Bolgnevllle '77 n 

2 

1 

2 



*** 

** 


Bolgnevllle ' 78 n 

6 

3 

2 



*** 

** 


Bolgnevllle '79 n 

3 

3 




*** 

** 


ISRAEL 

Kibbutz Boker '78° 

4 

1 



4 

** 

* 

** * 

SOOTH AFRICA 









.loodeplaat '78^ 

3 

1 



5 

** 

* 

* * 

AUSTRALIA 









Rutherglen '71 q 

1 

1 




* 

* 

*** 

Muzunbatenan '79 r 

1 

1 




* 

* 

*#* 

Murumba-enan '80 c 

1 

1 




* 

** 

a* 

Wagga Waf.ga '81 r 

1 

1 




* 

** 
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Table 1. Con't. 




Number of 


Measured 

Data in 

* 


Total 



Yield 




Data 

Sow Sow 


Phasic + Yld 

DM 

Soil 

Location 

Sets 

Var 1 Data Den 

Irr 2 

Devel Comp 

Partg 

Uatdir 


r.:XIC0 










Ciudad Obregon '73 s 

4 

1 



4 


* 



Ciudad Obzegon '74 c 

4 

1 



4 


* 

* 


U.S.A. 










North Dakota '79 c 

5 

5 




* 


** 


North Dakota '80 e 

4 

4 




* 

* 

** 


North Dakota '81 c 

4 

4 




* 

* 

** 


Sidney MN '77 v 

6 

3 


2 



* 

* 


Sidney MN ' 78 w 

6 

3 


2 



* 



Manhattan KS '82* 

2 

1 



2 

** 

** 

*** 


Hutchinson KS '80* 

1 

1 




** 

** 

** 


Hutchinson KS '81* 

1 

1 




** 

** 

* 


Carden City KS '81* 

2 

1 



2 

** 

** 

*** 


Garden City KS '82* 

4 

1 



4 

** 

** 

*** 


Bozeman MT '77 7 

i 

1 



i 

* 

★ 

* 


Lind HA '77* 

2 

1 


2 


** 

** 

* 


Pullman MA '73** 

2 

1 

2 



* 

* 

** 


Pullman HA '74** 

2 

1 

2 




** 

* 


Pendleton OR '77 bb 

2 

2 




** 

* 



Pendleton OR '80 bb 

8 

3 


3 


* 

** 

** 


Pendleton OR '81 bb 

1 

1 




* 

** 



Weston OR '81 bb 

3 

3 




* 

** 

* 


Temple TX '77 co 

21 

7 

3 




** 



Temple TX '83 dd 

8 

2 


4 


* * 


** 


Bus hi and TX '78** 

21 

2 

(2) 

2 

4 

#*■ 

** 


* 

Bushland TX '79** 

8 

1 



8 

*** 

** 

** 

** 

Bushland TX '80 ff 

1 

1 








Bushland TX '81 ff 

8 

1 



8 


** 

*** 


Bushland TX '82 £f 

8 

2 



4 


** 

*** 


Phoenix AZ '78** 

12 

2 



6 

* 

** 



Phoenix AZ '79“* 

19 

1 

5 


4 

** 

Hr 



Phoenix AZ '80 il 

17 

1 

5 


4 

** 

* 

**» 


Vernon TX '79 ff 

1 

1 




** 

* 

* 

* 


* Humber of asterisks reflects intensity of measurements . 

1 2 

Varieties in experiments! Irritation treatments. 

*Hatson et al. 1963 i b Pearman at al . 1978; c i d Baker 1979 i *Callather at al . 1916; 

f ; *Brooking 6 Kirby 1981; h Darwlnkel 1978; i Splertz & Ellen 1978; J Ellen t Spiertz 

1980; b Splertz l v.d.Haar 1978; ^Beinhauer, pars comm; ^ergermeier, pars comm; 

n Corre & Delacolle 1983; °Hochman 1979; p Meyer 1978; q Paltridge at al . 1972; 

c Armstrong, pars comm; *Sojka 1974; c Bauar 1980; v Aase 1978: w Black S> Aasa 1982; 

*Wagger 1983; 7 i * Johnson 1978; “Thill 1976; bb Klepper at al . 1983: ec Monk at 

al. 1979; dd 0teer; **Muslck & Dusak 1978: ££ Dusak, pars comm; bb ; il . 
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2. Independent data sets (153). These were run 
only on the. final version or, sometimes, on 
previous versions, but results of the run never 
affected the coding of the model. 

Methods 

This report is divided into (1) a report of the descrip- 
tive statistics to illustrate the general performance of 
CERES-Uheat and (2) a graphical comparison of model output 
with real world measurements of plant parts or growth factors 
over the growing season. The descriptive statistics consist 
of summary measures and difference measures. Criteria were 
selected by following the procedures recommended by Willmott 
(1982): 


observed mean (0) 
predicted mean (P) 

standard deviation of the observed values (S Q ) 

standard deviation of the predicted values 
(S p ) 

the intercept (a) and slope (b) of a simple 
regression of the form P - a + b * 0. The closer 
a to 0 and b to 1 the better the prediction 

the coefficient of determination of that regres- 
sion, describing to what extent the prediction 
accounted for variation of observed values (r^) 

Difference measures resulted partly from the above and 
included: 

mean bias error (MBE) 

- mean absolute error (MAE) 

- maximum absolute error (MAXAE) 

- minimum absolute error (MINAE) 

root mean square error (RMSE) 

All of these criteria are reported for the yield and 
yield components as well as for the dates of anthesis and 
physiological maturity, because matching of phenological 
stages is crucial for a valid run. 
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For the independent data sets, an "index of agreement" 
(d) was calculated according to the procedure of Willmot 
(1982). This calculation was only performed on yield data. 
All statistical computations were performed with the 
statistical package SAS (1982). 

The second and much more extensive part of the results 
are graphics, showing the proximity of the time course of 
model output and experimental measurements. Unfortunately, 
most data sets provided no information available on the 
variance of published measurements. This would be of certain 
interest because observed values are not true values in a 
statistical sense, but are mean values of a very small sample 
with more or less large intrinsic error. 


Results and Discussion 

Descriptive Statistics . Descriptive statistics based on 
the values shown in Tables 6 and 7 are reported separately 
for the two groups of data sets. Tables 2 and 3 are 
summarized measures for the dates of anthesis and physio- 
logical maturity, yield, grain weight, grain numbers per m^ 
(GPSM) , ear numbers per m^ , above ground dry matter, and leaf 
area index (LAI). Observed and predicted means are listed to 
show the difference between them. Standard deviations of 
observed and predicted values are listed to be checked for 
similar ranges. The intercept and slope of the regression 
equation with 0 as an independent variable provides informa- 
tion about the accuracy of predictions. Tables 4 and 5 
contain the difference measures for the same list of 
variables . 

Phasic Development . The accuracy of simulating the 
phasic development of a crop is crucial for the accurate 
simulation of crop growth and yield. Therefore, this part of 
the simulation received considerable attention. Figures 2-7 
show the scatter of predicted versus measured variables, 
including the 1:1 line +1 standard deviation and the matching 
regression. 

Regression equations for anthesis and physiological 
maturity show an excellent fit with slopes close to 1 and 
high regression coefficients. Observed means are consider- 
ably different for the two groups of data sets , the indepen- 
dent data sets reading anthesis and maturity 12 and 19 days, 
respectively, later. This reflects the fact that the origin 
of this subset of the data base is relatively far north, 
where the wheat crop season is extended into the latter part 
of the year. 


Difference measures in Tables 4 and 5 indicate small 
mean bias errors and mean absolute errors of five and seven 
days, respectively. The maximum errors for the two dates in 
the dependent subset occurred in the same run; therefore, the 
time available for grain filling was correct. When the date 
of anthesis was predicted too early, maturity came too early 
also. 


Yield. Grain Weight, and Grain Numbers/m ^. The indepen- 
dent data set is characterized by a higher mean yield caused 
by higher grain weight and slightly bigger grain numbers. 
This confirms that the data sets came from different 
locations. High yields were mostly from northern wheat 
growing areas. However, observed and predicted mean yields 
are very close resulting in small mean bias errors, llie mean 
absolute error is in the range of 22 percent of the observed 
mean yield and GPSM, but only 15 percent of the observed 
grain weight. 

Figures 2-4 demonstrate the scatter of predicted versus 
measured data points around the regression line, the 1:1 line 
and within the limits of +1 standard deviation of the 
observed mean. The percentage of data points outside of 
these limits is small as to the total number of observations . 
This demonstrates the model's reliability in simulating 
yield. Comparison of the regression to the 1:1 line proves 
the good fit, especially for the independent data sets. The 
regression coefficient, however, does not indicate a signifi- 
cant superiority of yield predictions over the rest of the 
data base (0.633 compared to 0.617). 

It is interesting to note that measured and predicted 
mean yields within a subset did not differ nearly as much as 
between the two subsets. 

Final Ear Number. Final Drv Matter, and Maximum LAI . 
Additional variables to be checked are final ear number, 
final dry matter, and maximum LAI observed during the growing 
season. Since the average experiment had about five measure- 
ments in the season, there is no guarantee that a measurement 
was taken at the peak of leaf area development. There are 
relatively few observations in total. Therefore, observed 
maximum LAI has to be considered with precaution. 

Despite these restrictions, the mean bias error is not 
very large. The mean absolute errors, however, amount to 40% 
and 95% of the observed mean. Also, the regression 
coefficient is not very high (Figure 7). (Please refer to 
the section where predicted and observed LAI are compared 
over the course of time for particular experiments to 
demonstrate the model's capability to simulate the canopy 
development. ) 
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Leaves comprise most of the total amount of above ground 
dry matter (Figure 6) . Since final dry matter is highly 
correlated to yield, it must be accurately simulated. The 
model overpredicts dry matter (Table 4) by about 380 g/m^, 
which is considerably less than 1 standard deviation (525.5 
g/m^). In the regression, the model accounts for 41.1 
percent of the variation in the real world. 

For the group of independent data, simulation of final 
ear number (Figure 5) resulted in a very small slope (0.1398) 
and corresponding coefficient of determination (r z - 0.014). 
The regression coefficient for the dependent data sets is 
0.40 and thus is higher than for dry matter. 

Difference Measures . While summary measures try tc 
describe the quality of simulation, difference measures try 
to locate and quantify errors. Difference measures for the 
same variables as above are listed in Tables 4 and 5. The 
mean bias error (MBE) is regularly conriderably smaller th 
the mean absolute error. The two values, being close, 
indicate the prediction to be biased in one direction, as is 
true for ear number and dry matter predictions with the 
independent sets . 

A negative MBE occurs when predictions are smaller in 
value than observations. Considerable errors of this kind 
are found with GPSM, ear number, and dry matter predictions 
of the dependent data sets, and with grain weight and date of 
anthesis of the independent sets. For all others, predic- 
tions are larger than observed values. 

Since the value of the MBE is related to the range of 
values, it is not reasonable to compare across variables. 
However, comparisons between the two groups of data sets 
might show the necessity of requiring independent data sets 
if it could be shown that the MBE and MAE of independent data 
were always larger. Unfortunately, while true for yield, it 
doesn't apply for all variables; maximum error in predicting 
yields is smaller in the independent set. 

Willmott (1982) indicated that the closer the root mean 
square error (RMSE) to 0, the better the model. Certainly 
this applies to the systematic error. The unsystematic RMSE 
should approach RMSE in the system: 

MSE - MSE + MSE 
s u 

The tendency of RMSE toward 0 is obvious, compared to RMSE. 
Since the mean is involved in the computation of this 
measure, comparisons across groups of data sets are not 
advisable. 
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The values of RMSE and RMSE, computed to predict yields 
for the dependent data set, are very similar: 25377 and 
25424, respectively. Thus, compared to RMSE of 1552, the 
model meets the above requirement. 

Willmott (1982) also suggests the computation of d, a 
quality measure for models, or index of agreemert: 


d - 1 


[ Z (P - 0.) / S 2 (|P'| + |0'|) ], 0 2 < d < 1 
i-1 i-1 


where P! - P. - 0 and 0! - 0. - 0 
i l li 


For a good model, d should approach 1. For CERES -Wheat yielu 
predictions, d equals 0.8825. This quality measure is mainly 
used to compare different models, but it provides additional 
useful information of model performance when used to evaluate 
a single model. 

Measured and Model Estimates of Crop Growth . To ensure 
accuracy, any simulation of plant parts should be checked 
intensively against real world observations. The appendix 
contains comparisons between model output with all available 
experimental measui ements . For some, real world data were 
obtained by digitizing graphs in publications, which adus 
uncertainty to the reliability of these data. These compari- 
sons should be viewed with a critical mind. Figures in 
Appendix A are in the order of locations . 

One problem with the data is that often only a few 
samples were taken during the growing season and for typical 
optimum curves like LAI or tiller number, the peak may be 
missing (see tiller number fo^ Rutherglen, Australia) . The 
experiment in Rutherglen is one of the few where root 
measurements were taken. It demonstrates good correspondence 
of compared data despite high variations in measurements and 
a slightly too high growth rate in the model run. Partition- 
ing was correct for this particular experiment as shown 
through a large number of samplings of root, stem, and ear 
weight. One discrepancy is that stem weight measurements 
seem to include leaves or sheaths , whereas the graph of the 
simulation shows stem weight only. 

Five, well-timed samples during the growing season 
demonstrated agreement with the tillering pattern at 
Murumbateman , Australia, in 1977. Biomass and LAI were 
simulated with high accuracy. The 1980 run at the same 
location did not perform quite as well, although the same 
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variety was grown. Data from Wagga Wagga, Australia, 1981, 
also showed excellent fit. 

The pattern of dry matter accumulation in response to 
different irrigation treatments was generally correctly 
simulated in Kibbutz Boker, Israel. The point of highest 
water stress is equally marked in the observations and the 
model simulations, although the differences between the two 
tend to increase toward the end of the growing period, 
causing the simulated LAI to increase too early and too much. 
Partitioning of assimilates to the ear, however, is very 
close to observations for those crops under late water 
stress . 

An irrigation experiment in Roodeplaat, RSA, shows 
considerable differences between the modelled and observed 
LAI. Plots receiving higher amounts of water (Irr. 3 and 5) 
produced much denser canopies which could be simulated well. 
In the dryland run (Irr. 1), the modelled LAI was too high, 
and it developed too late . 

In Ciudad Obregon, Mexico, another site where wheat is 
highly stressed by lack of water, the reverse effect was 
observed: LAI was generally underestimated, except for the 
plots experiencing the highest water stress where amount and 
timing were correct. 

In Bushland, Texas (1981 and 1982) , several different 
irrigation treatments for three varieties demonstrate 
excellent simulation of wheat growth and development with 
respect to dry matter, LAI, and tiller number. LAI shows the 
typical patcern for this area where wheat tillers vary 
greatly before winter due to warm fall temperatures, hence 
producing more coverage than further north where frosts 
usually kill many tillers and thus reduce leaf area. 
Therefore, these curves have two peaks. Some of the outliers 
in tiller number can be explained by the counting technioue 
used, which includes ever/ initiated tiller, even if not 
emerged yet. 

An experiment in Temple, Texas, in 1984 focused on the 
influence of sowing density on plant growth and tillering 
under restricted water conditions. The biomass simulation 
was excellent. LAI was sometimes larger than observed due to 
some problems in measuring leaf area before t^e leaves were 
fully dried. The model maintained a higher tiller number 
than observed in low sowing densities, but did not simulate 
the peak in the highest sowj density (640 plants/m^) . One 
hundred sixty plants/m z , a -v'mmon population for the area, 
gave excellent results . 
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A comparison with a similar experiment in a location 
with very different soil water availability in Lelystad, 
Netherlands, reveals interesting facts. The polder area of 
the Netherlands is characterized by almost unlimited water 
availability. In the model, the simulated crop produced an 
extremely high number of tillers, but were drastically 
reduced about 180 days after emergence. In reality, this 
decline in tiller number is in the same range, but it occurs 
earlier with higher populations. 

Other locations with the same climatic characteristics, 
like Flevoland and Uageningen, The Netherlands, demonstrate 
the excellent performance of CERES in cooler climates as 
well, though several runs suggested that it might perform 
better in the warm-to-subtropical climatic range. All of the 
biomass samples (5) taken at Flevoland in 1976, agree with 
the model's output. Biomass partitioning and LAI are 
represented by late season samples only, making profound 
comparison impossible, especially for LAI. Stem and ear 
weights agree reasonably well. Unusually high differences 
between simulated and measured leaf weights indicate a 
problem either with the definition of plant parts or units or 
something else. 

The same phenomena occurs in an experiment conducted the 
following year with the same variety, Lely, grown in 
Wageningen, The Netherlands. All the other details of this 
run show an excellent fit, as well as those with the variety, 
Maris Hobbit. With both, LAI simulation agreed perfectly 
with reality, including the timing of development, which is 
of particular interest, because Maris Hobbit was used in a 
series of experiments in France focusing on phenology. The 
French series was conducted for four years, wi".\ five 
varieties, at three locations' The growth details showed 
good fit with the model, v.hich is a profound test of the 
stability of the model's phenological predictions. This 
confirms how valuable experiment series are for model 
testing. 

Phenological observations complete with detailed 
measurements were taken in a three-year (1975-1977) series of 
experiments in Nottingham, England. Root weight measurements 
were considerably less than what the model had partitioned to 
the root system. Since root measurements are very difficult 
to take and therefore highly erratic, the model's functions, 
which work well for other locations, were not adjusted. 
Total biomass and its fractions gave excellent agreement with 
measurements. However, leaf weight was severely over- 
estimated in 1975, and in 1977 the crop developed too early 
and too fast in the model. 
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The less recent data sets are from Rothamsted, England, 
1960. The varieties Jufy and Atle are older, tail varieties; 
Capelle Desprey and Squarehead Master are new releases with a 
short culm. Therefore, only total dry matter and LAI could 
be compared, but with good results. The observed decrease in 
biomass at the end of the growing period does not appear in 
the simulations. The peak of LAI was generally under- 
estimated, except for Atle. 

The same tendency can be found in another highly 
productive crop in Ueihenstephan, Germany, in a record year: 
1983. For most of the growing season, coincidence between 
measured and modelled LAI prevailed, accompanied by excellent 
simulation of above-ground dry matter and tiller number. 

Several test runs were conducted with data from the main 
U.S. wheat producing areas: the Midwest and the Pacific 
Northwest. Experiments in Kansas (at Manhattan and at Garden 
City) focused on nitrogen fertilization; the irrigation 
aspect can also be checked with the non-nitrogen version of 
the model. In the Manhattan 1982 runs, the model recognized 
no water deficit for the non- irrigated crop. Therefore, 
model output was identical with the field. Measurements of 
details of the irrigated plots were slightly higher in the 
model, but they did not exceed estimated measurement 
variations. General agreement and timing is acceptable. 
Measured stem and leaf weights decreased more rapidly towards 
the end of the season than what the model allowed. Total 
measured biomass is larger, however, suggesting that a 
fraction, probably senesced leaves, is lacking. 

The 1981 Garden City experiment included two irrigation 
treatments; 1982 fovr. The results with the model are 
encouraging; dry matter was only slightly overe* ".imated, but 
LAI and tiller development were on target. There were hardly 
any differences between irrigation treatments two to four of 
the 1982 experiment in the model as well as in reality. In 
the tillering pattern there was no difference at all. 

The Sidney, Montana, series consists of two years of 
data for three varieties of Russian origin planted at two 
densities. The initial development of biomass was ahead in 
denser crops throughout winter and spring with the final 
results being identical. The model crop had 50% -80% more 
biomass, however, producing a much denser crop canopy than 
observed. The typical two-peak curve of leaf area develop- 
ment in the area with heavy winter-kill could not be checked, 
because observations started only in spring. Despite this 
lack of information, it is a very useful data set for testing 
the winter-kill and cold-hardiness routines of the CERES 
model. 
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The test runs for Lind, Washington, 1977, with two crop 
densities,- produced a very unusual pattern of leaf area 
development. Despite double the initial population, the 
course of tiller development in the model was almost 
identical to the low density crop where final tiller number 
coincided with measurements. The high density crop produced 
a considerably higher number of tillers than simulated. 
Since few late measurements were taken, validation of the LAI 
component is difficult. It seems the model was overesti- 
mating LAI at this location, as winter temperatures killed 
most of the leaf cover. 

Similar results, although less pronounced, were found 
for Pullman, Washington. With the 1973 data, CERES' 
simulations coincided perfectly with observations for early 
and late sowing; for 1974 the simulation was good only for 
the late sowing data. There is some confusion, though, with 
the data, as observed LAI and dry matter of the early sowing 
were unusually high. 

The Experiment Station at Pendleton, Oregon, furnished 
valuable data from diverse experiments. Detailed measure- 
ments from 1977 allowed us to check the correct simulation 
for different plant parts. For unknown reason, however, ear 
weight was far below the field- observations . Leaf area was 
simulated to be twice as high as reality. The tiller number 
comparison is impressive, although it is not certain that 
plotted tiller counts represent the peak of tillering. The 
final increase in tiller number is a very unusual phenomena. 

Biomass weights from a 1980 experiment in Pendleton and 
Weston (Oregon) with two sowing dates and three sowing 
densities, though few, confirm the accuracy of simulation. 
This was reinforced by the 1981 Pendleton run. The 
differences between the modelled and measured values are 
acceptable . 

The CERES model, developed from a sail water balance 
subroutine, had been tested and validated. Thus, few tests 
were done with soil water measurements. 
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Figure 2. 


Predicted Versus Measured Yield. 



Figure 3. Predicted Versus Measured Grain Weight. 
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Predicted Car Number pur sqm Predicted Grain Number per sqm 
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Figure 5. Predicted Versus Measured Number of Ears. 



17 
















Table 2. Summary Measures for Independent Data Sets. 


Variable 

Unit 

n 

0 

p 

S o 

S P 

P - a 
a 

+ b*0 
b 

Re t 

Coeff 

Anchasls 

days 

82 

1*4.6 

1*4.2 

30.3 

28.4 

14.6 

0.8960 

. 912 

Phys Mat 

days 

76 

176.0 

176.6 

28.4 

28.5 

4 . 7 

0 . 9763 

. 947 

Yield 

kg /ha 

157 

55*7.0 

56*6.8 

2219.0 

2578.2 

582.2 

0.9166 

.633 

Crain Uc 

me /kerne 1 

lit 

38.3 

36.4 

9.3 

9.4 

13.4 

0.6234 

. 448 

CPSM 

* 

138 

1*594.6 

15161.2 

5115.2 

5811.9 

5507.3 

0.6715 

.376 

Ear * 

i 

122 

524.0 

6*8.7 

152.9 

190.2 

600 . 3 

0 .1398 

.014 

Dry Matter g/m 2 

76 

1179.2 

1559.1 

525.5 

412.2 

791 . 8 

0.5052 

. *11 

LAI 

-- 

5* 

4 . 4 

5.2 

2.6 

1.6 

3 . 3 

0.4342 

. 260 


Table 3. 

Summary Measures for Dependent 

Data Se t s . 




Variable 

Unit 

n 

Meas 

Mean 

Pred 

Mean 

s « 

S P 

P - a 
a 

+ b*0 
b 

Reg 

Coeff 

Anthesls 

days 

113 

132.9 

135.7 

62.1 

62.8 

2.46 

1 . 0021 

. 983* 

Phys Mat 

days 

111 

157.2 

151.9 

50.3 

54 . 0 

1.56 

0.9735 

. 9*16 

Yield 

kt/ha 

130 

4359.6 

4376.0 

1892.7 

1892.9 

950.7 

0 . 7857 

. 6171 

Grain Ut 

mg/ka rnel 

91 

31.9 

33.1 

6.2 

8.4 

8.3 

0 . 7598 

. 3*67 

CPSM 

t 

91 

13837.9 

1301 • .8 

4144.3 

4*77.1 

5705.8 

0.5596 

. 3167 

Ears 

f 

77 

558.2 

527 . 3 

192.2 

270.1 

185.3 

0.7818 

. 3998 

2 

Dry Matter g/m 

85 

1396.4 

1359.6 

529.9 

418.4 

1033.1 

0.3296 

. 33.’ 5 

LAI 

-- 

37 

4.0 

4.9 

1.9 

2.4 

8.26 

-0.2957 

.0395 


Table 4. 

Difference 1 

Measures 

for Independent Data 

Sets . 








Max 

Min 


Variable 

Unit 

n 

MBE 

MAE 

Error 

Error 

RMSE U 

Anthesls 

days 

82 

-0.4 

5.67 

32.0 

0 . 

8.5 

Phys Mat 

days 

76 

0.6 

5.1 

18.0 

0 . 

6.6 

Yield 

kg /ha 

157 

99.8 

1272.8 

3468.0 

6 . 

1552.0 

Grain Wt 

mg /kernel 

144 

-1.9 

5.8 

28.4 

0 . 

6.5 

CPSM 

#/m 2 

138 

566.6 

3*44.1 

17*46.0 

17. 

4437 . 5 

Ear # 

#/m 2 

122 

124.7 

196.0 

938.0 

1 . 

179.6 

Dry Matter 

g/m 2 

76 

379.9 

352.0 

1324.0 

2. 

320 . 0 

LAI 

-- 

54 

0.8 

1.8 

5.7 

0 . 

1 . 9 


Table 5. 

Difference 

Measures 

for Dependent Data 

Sets . 



Variable 

Unit 

n 

MBE 

MAE 

Max 

Error 

Min 

Error 

RMSE U 

Anthesls 

days 

113 

2.8 

6.8 

24.0 

0 . 

8.1 

Phys Mat 

days 

111 

-5.4 

7.7 

24.0 

0 . 

12.2 

Yield 

kg /ha 

_.»J 

16.4 

948.3 

4066.0 

19. 

1175.8 

Crain Wt 

mg /kernel 

91 

1.2 

5.1 

22.2 

0.1 

6.5 

CPSM 

f/m 2 

91 

-821.1 

2907.1 

11059.0 

72. 

3425.6 

Ears 

#/m 2 

77 

-30.9 

150.4 

563.0 

1 . 

185.3 

Dry Matter 

g/o 2 

85 

-36.8 

354.0 

1314.0 

1 . 

247.8 

LAI 

-- 

37 

0.9 

3.8 

9.6 

0 . 

2.8 


Table 6. Listing of Final Results of Independent Data Sets. 


YIELD 

GRAIN 

WT. 

GRAINS /SQM 

EARS /SQM 

MAX. 

LAI 

DRY 

MATTER 

0 

P 

0 

P 

0 

P 

0 

p 

0 

P 

0 

P 

1700. 

1319. 

23.3 

20.5 

7205. 

6*30 . 

275. 

*22- 

1 . 1 

2.0 

*48 . 

*86 . 

3575. 

3600 . 

32.0 

34.2 

11281 . 

10521 . 

357 . 

365. 

4.0 

2.1 

775 . 

826 . 

2110. 

299* . 

26.5 

31.6 

7945. 

9*79. 

373. 

3*8. 

2.0 

1.9 

5*3 . 

752. 

1826. 

2803. 

29.0 

35.4 

6169. 

791*. 

298 

4*3 . 

0.0 

2.8 

409. 

850 . 

1180. 

3068. 

20.7 

23.5 

5636. 

13056. 

0. 

*63. 

1 . 0 

2.6 

275. 

8*5. 

2933. 

3550 . 

29.6 

28.1 

99*6. 

12623. 

0 . 

513. 

3.0 

2.8 

662 . 

891 . 

256*. 

3879. 

28.7 

30.0 

8995. 

12935. 

0. 

*61. 

1 . 9 

2.8 

523 . 

935. 

2*8*. 

3*6*. 

31.1 

34.7 

8017 . 

9989. 

0. 

*02. 

3.5 

2.4 

635. 

8*5. 

2239. 

3799. 

24.8 

29.3 

91*5. 

129*5. 

0 . 

517. 

2.3 

2.9 

595. 

1010 . 

3*2* . 

3578. 

26.5 

29.3 

129*5 . 

12225. 

0 . 

*73 . 

3 . 6 

2.6 

999. 

9*5. 

2760. 

*031 . 

23.6 

31.7 

11726. 

12710. 

0 . 

530. 

2.8 

3 . 4 

611 . 

1053 . 

2330 . 

3111. 

26.8 

30.9 

3689 

10052. 

0. 

*1* . 

3.7 

2.8 

631 . 

908. 

6960. 

8762. 

0.0 

*3.3 

0. 

20219. 

518. 

604 . 

0.0 

5.5 

0 . 

1892. 

5*70. 

7679. 

0.0 

*4.7 

0 . 

17165. 

*80 . 

575. 

0.0 

5.2 

0 . 

1690. 

65*0. 

8755. 

0.0 

*1.8 

0 . 

209*5 . 

0 . 

567. 

0.0 

5.6 

0 . 

1945. 

6340. 

8*0* . 

*8.6 

*1 . 7 

0 . 

20138. 

0 . 

609. 

0.0 

5.2 

0. 

1920 . 

5620. 

6660. 

48.0 

39.6 

0. 

16788. 

*85 . 

55* . 

0.0 

4.7 

0 . 

1505. 

7080. 

7369. 

51.0 

40.7 

0. 

1809*. 

*25. 

585. 

0.0 

* . 7 

0 . 

1568 . 

6990. 

7966. 

0.0 

** . 2 

0 . 

18023. 

623. 

530. 

0 . 0 

4.5 

0 . 

1678. 

5810. 

805* . 

38.4 

39.8 

0 . 

20227. 

47* . 

56*. 

0.0 

4.8 

0. 

1756. 

6500. 

7872. 

36.1 

37.3 

0 . 

21128. 

6*6 . 

603. 

0.0 

5.4 

0 . 

1799. 

8300. 

8969. 

53.7 

*3.7 

0 . 

20533. 

467. 

530. 

0 . 0 

5.4 

0 . 

1932. 

7810. 

8042. 

*9.0 

37.2 

16138. 

21597. 

492- 

642. 

0.0 

5.6 

0 . 

1980 . 

8690 . 

7132. 

4* . 3 

35.9 

19918. 

19882. 

57* . 

575. 

0.0 

5.0 

0 . 

1743 . 

10760. 

8203 . 

*8.4 

*1.3 

22784. 

19862. 

6*0 

655. 

0.0 

6.1 

0 . 

1988 . 

10900. 

8930. 

52.2 

** . 3 

2116* . 

20157. 

67*. 

539. 

0.0 

5.* 

0 . 

1904 . 

9560. 

7600. 

5* . 3 

*6.5 

17970. 

1634*. 

599. 

*50. 

0 . 0 

3.8 

0 . 

1526. 

95*0. 

7836. 

*2.1 

A* . 0 

2300* . 

17818. 

695. 

595. 

0.0 

5.7 

0 . 

1791. 

9510. 

714* . 

*8-5 

** . 7 

19880 . 

1599* . 

568 . 

520. 

0 . 0 

4 . 7 

0. 

16*2. 

9960. 

7*55. 

*9.1 

*2 . 9 

20376. 

1738* . 

532. 

508. 

0.0 

* . * 

0 

1*79. 

9280. 

7879. 

*0.8 

*1.3 

22616. 

19079. 

57* . 

525. 

0 . 0 

5.6 

0 , 

1856 . 

3361. 

2*03. 

.0 

20.0 

0 . 

12017. 

26* . 

507. 

5 

2.8 

773. 

922. 

*213. 

1*76. 

.0 

20.0 

0 . 

7381 . 

*8* . 

626 . 

3.2 

6.1 

1576, 

1197. 

*47*. 

1*27. 

.0 

20.0 

0 . 

7136. 

515. 

6*5. 

4.2 

9.1 

1497 , 

1478. 

*712. 

12*4. 

.0 

20.0 

0. 

6222. 

711 . 

659. 

9.6 

10.2 

2518 

15*1 . 

0 . 

1170 . 

.0 

20.0 

0 . 

5851. 

37* . 

57* . 

2 . 1 

3 . 3 

955 

860 . 

3102. 

659. 

.0 

20.0 

0 . 

3296. 

560 . 

664 . 

* . 0 

6.9 

1178. 

1151 . 

7391 . 

637. 

.0 

20.0 

0 . 

3186. 

960 . 

670. 

4.5 

10.1 

1505 

1*34 . 

3972. 

595. 

.0 

20.0 

0 . 

2975. 

*95. 

660 . 

12.2 

11.3 

1922 

1*9* 

3116. 

17*5. 

32.1 

37 . 8 

9707. 

4620. 

0. 

593. 

3.3 

3 . 1 

1008 . 

913 . 

4561 . 

A 822. 

32.3 

37.8 

1*121 . 

12768 . 

0. 

620. 

3.6 

3.1 

1285 

1*06. 

2633. 

1338. 

31.0 

33 . 3 

8*9* . 

4620. 

0. 

593 . 

3 . 3 

3.1 

95 7 

768. 

5945. 

*322. 

36.9 

37.8 

16111 . 

12768. 

0. 

620. 

Z . f 

3.1 

1184 

1406 . 

3102. 

3181. 

26.9 

24.5 

11532. 

12979. 

0 

802. 

2.6 

3.5 

1217 

1203 

*1*5. 

*853. 

21.3 

23.6 

19*60. 

20576. 

0. 

902. 

3.9 

6.‘ 

1303 . 

1723 . 

2277. 

110* . 

23.7 

22.0 

9608. 

5028 . 

0. 

794. 

2.1 

3.5 

687. 

713 . 


20 
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Table 6. Continued. 


YIELD 

GRAIN 

WT . 

GRAINS/ SQM 

EARS /SQM 

MAX. 

LAI 

DRY 

MATTER 

0 

P 

0 

p 

0 

P 

0 

P 

0 

P 

0 

P 

5833 . 

5043 . 

25.6 

24.5 

22785. 

20576. 

0 . 

902. 

3 . 9 

6.6 

1771. 

1881 . 

2*77 . 

704. 

28.4 

35.4 

8721. 

1990. 

0. 

726. 

2.8 

5.4 

606. 

751 . 

3752. 

4095. 

25.5 

36.8 

14713. 

11119. 

0. 

826. 

5.3 

5.5 

838. 

13*2 . 

*12* . 

2778. 

31.8 

41 . 1 

12968. 

6761. 

0 . 

726. 

5.3 

5 . 4 

762. 

1136 . 

4009. 

818. 

29.3 

41 . 1 

13682 . 

1990 . 

0 . 

726. 

5.3 

5.4 

700. 

884 . 

4550. 

5456. 

31 . 4 

41 . 1 

14490 . 

13277 . 

0 . 

826. 

5.3 

5.5 

796. 

1569. 

4337 . 

4569. 

30.6 

41.1 

14173 . 

11119. 

0 . 

826. 

5.3 

5.5 

838. 

1426. 

3991. 

2778. 

31.1 

41.1 

12832. 

6761 . 

0. 

726. 

5.3 

5.4 

706. 

1139. 

5689. 

5456. 

34.7 

41.1 

16394. 

13277. 

0 . 

826. 

3 . 9 

5.5 

7*6. 

1569. 

8967. 

7662. 

47.1 

58.1 

19042. 

13197 . 

422. 

540 . 

. 0 

4.6 

0 . 

1772. 

8351 . 

6406. 

43.2 

45.9 

19350 . 

13951 . 

386. 

425. 

.0 

4.6 

0 . 

1740 . 

7694 . 

7333. 

25.5 

39.9 

30197 . 

18381. 

466. 

561. 

.0 

4.6 

0 . 

1766. 

4950. 

4267. 

38.0 

41.7 

11720. 

10228. 

0 . 

267. 

. 0 

3.8 

0 . 

13*5. 

2920. 

3963. 

35.0 

43.0 

7390. 

9226. 

0. 

279. 

.0 

3.7 

0 . 

1255. 

3720. 

4267 . 

33.0 

41.7 

10180 . 

10228. 

0 . 

267 . 

.0 

3.8 

0 . 

1292. 

3830 . 

4219. 

32.0 

38.0 

10760 . 

11110. 

0. 

277. 

. 0 

3 . 4 

0 . 

12*3 . 

5650. 

2461. 

53.3 

24.9 

10600. 

9871 . 

0 . 

518. 

0.0 

5.5 

1212. 

1298. 

6640 . 

3761. 

47.6 

26.2 

13950. 

14356. 

0 . 

505. 

0.0 

5.8 

1198 . 

1368. 

5500. 

7899. 

43.0 

54.4 

12825. 

14514 . 

*75. 

546. 

10.2 

7.1 

1233 . 

2121 . 

6100. 

9232. 

44 . 0 

39.6 

13775. 

233*2. 

*75. 

798. 

11.6 

5.9 

1195. 

2173 . 

5900 . 

8231. 

40.0 

45.5 

14520. 

18102. 

440 . 

903 . 

8.2 

5.9 

1151 . 

2100. 

4062. 

2428. 

42.0 

23.0 

9672. 

10578. 

372. 

5*5. 

3.7 

8.3 

851 . 

1628. 

4072. 

4009. 

42.3 

20.7 

9628. 

19395. 

332. 

779. 

3.4 

7.2 

802 . 

1720. 

4250. 

3750. 

36.8 

24.9 

11550. 

150*1. 

350. 

523. 

3.2 

7.2 

814 . 

1711. 

6794 . 

7991. 

33.0 

40.6 

20587. 

19686. 

67J. 

606. 

0.0 

6.1 

2158. 

1958. 

6922. 

8599. 

40.5 

38.5 

17091. 

22330. 

79*. 

720. 

0.0 

6.1 

1/ 40 . 

1994 . 

7050. 

7834. 

43.6 

40 . 6 

16170. 

19292. 

652. 

694. 

0.0 

6.1 

183*. 

1948 . 

4865. 

7783. 

28.3 

39.3 

15106. 

19811 . 

506. 

609. 

11.8 

6.2 

1*15. 

1951 . 

6000. 

7996. 

32.9 

40.9 

18215 . 

19539. 

600. 

720 . 

0.0 

6.1 

1602. 

1950. 

7886. 

5886. 

35.6 

29.3 

22102. 

20070. 

682. 

661. 

0.0 

4.7 

1891 . 

1*28. 

7088. 

5722. 

38.8 

32.5 

18270 . 

17605. 

599. 

635. 

0.0 

4 . 7 

1739. 

1*05. 

6930. 

5722. 

39.2 

32.5 

17661 . 

17605. 

5*6. 

635. 

0.0 

4 . 7 

1*60. 

1*05. 

7991 . 

5890. 

34.8 

28.8 

22995. 

"C422 . 

640. 

635. 

0.0 

4 . 7 

1772 

1428. 

6570 . 

8067. 

46.0 

40.1 

14900 

20098. 

429. 

387. 

4.3 

5.1 

1682 

1789. 

7590 . 

6835. 

50.8 

51 . 9 

164C 

13167 . 

462 . 

*50. 

5.2 

5.1 

1683 

1685. 

8580. 

6619. 

42.1 

27.2 

20552. 

243*0. 

538. 

390. 

0.0 

5.4 

1787 

1737 . 

8260. 

9731. 

41.4 

48.7 

19970 . 

20000. 

635. 

775. 

0.0 

5.7 

1810 

2020. 

9390 . 

9857. 

56.3 

47.2 

16775 . 

20880 . 

510. 

591. 

5.0 

5.5 

1700 

2088. 

4480 . 

7 535. 

52.3 

58.1 

8619. 

13140 . 

650. 

530. 

3 . 3 

5.0 

940 

1516. 

6310 . 

4492. 

53 . 7 

42.3 

11761 . 

10611. 

420. 

574 . 

0.0 

5.5 

980 

1417 . 

74' 0 . 

8147. 

49.2 

39.0 

15229. 

20880. 

490. 

591. 

0.0 

5.5 

1250 

1901 . 

7090. 

8702. 

38.8 

35.1 

18277 . 

2*773. 

518. 

981 . 

0.0 

7.8 

0 

2320. 

5560. 

7963. 

40.3 

34.4 

13810. 

23163. 

440. 

9*6. 

0.0 

7 . * 

0 

2037 . 

6810. 

8661 . 

49.7 

43.9 

13696. 

19716. 

391 . 

1308. 

0.0 

7.5 

0 

2298. 

5710. 

7925. 

50.5 

44 . 6 

11313. 

17762. 

288. 

1226. 

0.0 

7.4 

0 

2038. 

6440. 

8008. 

33.6 

38.4 

19163. 

20878. 

645. 

808. 

0.0 

«. 7 

0 

2042. 

6930. 

7758. 

41.8 

.9.0 

17431. 

15844. 

635. 

792. 

0.0 

7.0 

0 

2009. 
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T»bl« 6 . Continued. 


YIELD 

GRAIN 

WT. 

GRAINS/SQM 

EARS /SQM 

MAX. 

LAI 

DRY 

matter 

0 

P 

0 

P 

0 

P 

0 

P 

0 

p 

0 

p 

5920 . 

8049 . 

335 

37 . 9 

17668 . 

21242 . 

568 . 

786 . 

0.0 

6.0 

0 . 

1953 . 

5180 . 

7200 . 

40.8 

44.2 

12669 . 

16284 . 

543 . 

782 . 

0.0 

5.9 

0 . 

1874 . 

8740 . 

8594 . 

35.1 

33.9 

24979 . 

25341 . 

760 . 

838 . 

0.0 

6 . 7 

0 . 

2159 . 

7720 . 

7995 . 

34.6 

31.8 

22304 . 

25158 . 

562 . 

828 . 

0.0 

5.9 

0 . 

2020 . 

6460 . 

5027 . 

28.8 

20.0 

22405 . 

25111 . 

737 . 

838 . 

0.0 

6. 7 

0 . 

1757 . 

5010 . 

4474 . 

24.0 

20.0 

20905 . 

22371 . 

534 . 

828 . 

0.0 

5.9 

0 . 

1626 . 

7910 . 

8671 . 

46.4 

44 . 1 

17039 . 

19667 . 

629 . 

1071 . 

0.0 

6.3 

0 . 

2117 . 

6850 . 

5190 . 

42.2 

26.4 

16241 . 

19665 . 

589 . 

1071 . 

0.0 

6.3 

0 . 

1758 . 

5191 . 

8448 . 

39.2 

38.4 

13100 . 

21999 . 

413 . 

858 . 

0.0 

6.0 

0 . 

2082 . 

5590 . 

8700 . 

38.1 

39.5 

14502 . 

21999 . 

433 . 

858 . 

0.0 

6.0 

0 . 

2107 . 

8890 . 

9021 . 

40.0 

37.4 

22225 . 

24093 . 

759 . 

870 . 

0.0 

6.6 

0 . 

2190 . 

8290 . 

8935 . 

52.1 

48.7 

15911 . 

18362 . 

499 . 

1143 . 

0.0 

6.3 

0 . 

2148 . 

8630 . 

8418 . 

48.5 

46.6 

17794 . 

18050 . 

603 . 

868 . 

0.0 

6.9 

0 . 

2143 . 

8530 . 

8221 . 

47.0 

41.9 

18149 . 

19630 . 

626 . 

1062 . 

0.0 

6.6 

0 . 

2111 . 

7546 . 

7559 . 

43.2 

35.2 

17467 . 

21452 . 

683 . 

918 - 

0.0 

6.7 

0 

2017 . 

6992 . 

6804 . 

51.1 

43.4 

13677 . 

15687 . 

614 . 

909 . 

0.0 

6.8 

0 

1968 . 

4810 . 

6632 . 

33.3 

28.6 

10589 . 

23223 . 

460 . 

849 . 

0.0 

6.1 

0 

1843 . 

4550 . 

5667 . 

34.2 

26.0 

9792 . 

21760 . 

401 . 

790 . 

0.0 

5.1 

0 . 

1503 . 

4420 . 

6725 . 

41.1 

37.7 

8423 . 

17862 . 

300 . 

1115 . 

0.0 

5 . 9 

0 . 

1833 . 

4010 . 

5697 . 

41.7 

34 . 6 

6715 . 

16453 . 

288 . 

997 . 

0.0 

4.8 

0 . 

1458 . 

5600 . 

8381 . 

36.5 

36.3 

15263 . 

23117 . 

478 . 

766 . 

0.0 

5.4 

0 . 

1877 . 

5700 . 

7463 . 

36.0 

35.4 

15814 . 

21108 . 

548 . 

772 . 

0.0 

5.1 

0 

1724 . 

6780 . 

8795 . 

45.1 

46. 4 

15013 . 

18961 . 

391 . 

655 . 

0.0 

6.6 

0 . 

2066 . 

6080 . 

8665 . 

45.5 

46.3 

13367 . 

18701 . 

377 . 

628 . 

0.0 

5.1 

0 

1770 . 

7710 . 

9064 . 

38.1 

42.7 

19922 . 

21208 . 

686 . 

755 . 

0.0 

6.6 

0 

2079 . 

7830 . 

8657 . 

38.1 

42.6 

18976 . 

20311 . 

606 . 

718 . 

0.0 

5.1 

0 

1794 . 

7090 . 

8495 . 

45-6 

52.2 

13553 . 

16260 . 

513 . 

740 . 

0.0 

6 . / 

0 

2048 . 

6480 . 

8148 . 

40.5 

50.8 

15980 . 

16034 . 

664 . 

708 . 

0.0 

5 . 1 

0 

1744 . 

7360 . 

7999 . 

42.6 

44 3 

17446 . 

18067 . 

494 . 

628 . 

0.0 

5.2 

0 

1674 . 

7010 . 

8058 . 

43.5 

41.2 

16168 . 

19571 . 

463 . 

686 . 

0.0 

4 . 7 

0 

1599 . 

7280 . 

7652 . 

47.0 

49.5 

15479 . 

15459 . 

484 . 

701 . 

0.0 

5.1 

0 

1626 . 

4910 . 

4356 . 

39.0 

36.7 

11190 . 

11880 . 

0 . 

326 . 

0.0 

4 . 4 

0 

1440 . 

4130 . 

4706 . 

33.0 

40.8 

11350 . 

11535 . 

0 . 

286 . 

0.0 

3 . 5 

0 

1363 . 

3730 . 

4356 . 

33.0 

36.7 

10210 . 

11880 . 

0 . 

326 . 

0.0 

4 . 4 

0 

1409 . 

2710 . 

3574 . 

32.0 

32.7 

7570 . 

10946 . 

0 . 

295 . 

0.0 

4 . 3 

0 

1262 . 

1950 . 

2386 . 

39.1 

34.9 

4987 . 

6839 . 

i 48 . 

638 . 

0.0 

3 . 5 

0 

887 

1900 . 

2042 . 

37 . 9 

34 . 9 

5015 . 

5854 . 

170 . 

644 . 

0.0 

3.6 

0 

892 . 

2372 . 

2366 . 

32.1 

27.4 

7390 . 

8637 . 

263 . 

639 . 

0.0 

3.5 

0 

892 

2217 . 

2234 . 

31.3 

29.7 

7084 . 

7518 . 

308 . 

643 . 

0 . 0 

3.6 

0 

896 

2450 . 

2281 . 

28.3 

29.4 

8641 . 

7766 . 

370 . 

642 . 

0.0 

3.4 

0 

890 

2352 . 

1885 . 

27.9 

29.4 

8433 . 

6417 . 

499 . 

640 . 

0.0 

3.8 

0 

889 

3455 . 

4304 . 

32.1 

35.6 

10763 . 

12088 . 

379 . 

655 . 

0.0 

4.0 

829 

1414 

3460 . 

4378 . 

33.3 

35.6 

10651 . 

12297 . 

433 . 

664 . 

0.0 

4 . 7 

842 

1520 

4277 . 

4848 . 

32.4 

31.4 

13200 . 

15451 . 

440 . 

656 . 

0.0 

4 . 0 

1044 

1439 

4147 . 

4953 . 

31.4 

31.4 

13209 . 

15786 . 

518 . 

667 . 

0.0 

4 . 9 

1024 

1566 

3888 . 

4228 . 

30.3 

30.0 

12831 . 

14109 . 

578 . 

656 . 

0.0 

4.0 

1022 

1421 

3520 . 

4320 . 

29.7 

30.0 

31862 . 

14416 . 

674 . 

667 . 

0.0 

4.9 

986 

1553 
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Table 6. Continued. 


YIELD 

GRAIN 

UT. 

GRAINS /SQM 

EARS / SQM 

MAX. 

LAI 

DRY 

MATTER 

0 

P 

0 

P 

0 

P 

0 

P 

0 

P 

0 

P 

4830 . 

2772. 

40.1 

28.0 

12125. 

9900. 

548. 

452. 

0.0 

2.2 

883 . 

693 . 

2172. 

1278. 

32.0 

20.0 

10353. 

6388. 

1094. 

740 . 

2.0 

7 . 2 

1279. 

1309. 

2032. 

1316. 

33.3 

20.0 

8008. 

6579. 

654. 

742. 

3.4 

7.5 

931. 

1305. 

3738. 

3763. 

0.0 

20.0 

0 . 

18813. 

621. 

555. 

5 . 5 

5.5 

1300 . 

1537. 

4642. 

4597 . 

0.0 

20.0 

0. 

22985. 

625. 

509 . 

4 . 1 

4.8 

1400 . 

1426. 

8778. 

3339. 

35.4 

33.1 

24520. 

25161. 

757. 

575. 

87 

5.0 

3389. 

2065. 

7109 

8063. 

33.6 

35.3 

2103L. 

22866. 

589. 

504 . 

4.9 

4 . 4 

1778. 

1766 . 

2420 . 

1070 . 

13 . 4 

20.0 

21270. 

5351. 

600. 

814 . 

3.0 

6 . 6 

1361 . 

1228. 

2250 , 

1087. 

12.5 

20 . 0 

17900. 

5436. 

800. 

814 . 

2.7 

6 . 6 

1215. 

1227 . 

7649. 

6315. 

55 . 9 

49.9 

13683. 

12662. 

449. 

508. 

0.0 

4 . 0 

0. 

1605 . 

6104 . 

6483 . 

54 . 1 

49.9 

11277. 

12999. 

467. 

542. 

0.0 

5.1 

0 . 

1761 . 

5897 . 

6523. 

55.4 

49 . 9 

10648. 

13078. 

488. 

545. 

0.0 

5 . 5 

0 . 

1814 

6900 . 

6639. 

55.3 

53.1 

12485. 

12502. 

357. 

459. 

0 . 0 

3 . 7 

0 . 

146? . 

6259 

6414. 

53.4 

49.8 

11718. 

12879 . 

413. 

518. 

0.0 

4 . 4 

0 

1545. 

5992 

6265. 

55.0 

48.1 

10905. 

13024. 

454 . 

541 . 

0 .U 

4 . 7 

0 

1596. 

0466 

3912. 

45.9 

43 . 4 

14102. 

13615. 

362. 

411 . 

0.0 

4 . 4 

0 

1530. 

6472 

6564. 

34.6 

36.6 

18835. 

17939. 

437 . 

537. 

0.0 

4 . 4 

0 

1544. 

9665 

6485. 

51.9 

54 . 1 

18622. 

11993. 

632. 

536. 

0.0 

5.0 

0 

1781. 
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Tabla 7. Listing of F'LnaL Rosales of Dopondsnt Data Sats. 


YIELD 

GRAIN 

WT. 

GRAINS/ SQM 

EARS /SQM 

MAX. 

LAI 

DRY 

MATTER 

0 

P 

O 

P 

0 

P 

0 

P 

0 

P 

O 

P 

4825. 

7843. 

45.4 

64 . 5 

10670 . 

12164 . 

485. 

465. 

7.7 

6.1 

1500 . 

1914 . 

2320. 

3002. 

37.0 

40.3 

5867 . 

7442 . 

402. 

546. 

7 . 9 

7 . 1 

1250. 

1847. 

4070. 

4800. 

43.2 

45.5 

9325. 

10543 . 

374. 

400. 

5.1 

4.0 

1135. 

1382. 

2710. 

3585. 

32.3 

36.4 

8175 . 

9859. 

415. 

499. 

4.8 

5 . 0 

1050 . 

1521 . 

6450. 

6469. 

0.0 

32.2 

0 . 

20064. 

0 . 

701. 

4.0 

5.3 

1252. 

1725. 

3400 . 

3594. 

0.0 

20.8 

0. 

17263. 

0. 

689. 

0.0 

6 . 5 

0. 

1543. 

5030. 

7961. 

0.0 

39.1 

0. 

20379. 

413. 

733. 

0.0 

5.3 

0. 

1909 

2820. 

5058. 

.0 . 3 

38.3 

6997 . 

13220. 

118. 

155. 

0.0 

1.6 

707 . 

884 . 

6440. 

7065. 

41 . 7 

39.2 

15443. 

18020. 

272. 

264 . 

0.0 

3 . 2 

1531 . 

1331 

7200 . 

7574 . 

44 . 6 

39.9 

16143 . 

18996. 

322 . 

311 

0.0 

4 . 0 

1647 . 

1480 . 

8430. 

7893. 

44 . 1 

40 . 4 

19115 . 

19534 . 

430 . 

359. 

0.0 

4 . 7 

1968. 

1616 . 

8340. 

8065. 

42.3 

40.8 

19716. 

19788. 

490 . 

383. 

0 . 0 

5.2 

2024 . 

1753 . 

8480. 

8078. 

41.0 

40.8 

20683 . 

19809 . 

582 . 

400 . 

0.0 

5 . 2 

2181 . 

1814 . 

8900 . 

8180 . 

39.0 

41.1 

22820 . 

19898. 

777. 

800 . 

0.0 

5.7 

2376. 

1870 . 

4440. 

3125. 

38 . 2 

25.6 

11642 . 

12227 . 

433 . 

564 . 

4.0 

4 . 8 

109-. . 

1196. 

2580 . 

2349. 

35.5 

25.9 

7319 . 

9060. 

308. 

587 . 

2.2 

5.0 

717 . 

1107 . 

3350 . 

1752. 

21.3 

20.0 

15830. 

8761. 

581 . 

577. 

4 . 9 

4 . 9 

1115. 

1093 . 

4344 . 

4707. 

34.0 

39.1 

12800 . 

12053 . 

466. 

374. 

6 . 5 

5 . 0 

1492. 

1545. 

8711 . 

6876. 

44 . 0 

43 . 2 

19934 . 

15906 . 

511 . 

728. 

0.0 

5 . 6 

0. 

1692. 

3537 . 

1992. 

24.0 

20.0 

14738. 

9960 . 

391 

714. 

0.0 

4 . 8 

0 . 

1096 . 

2508. 

1368. 

20.6 

20.0 

12175. 

6340. 

372. 

680. 

0.0 

5.2 

0 , 

980 . 

3009. 

1654. 

23.0 

20.0 

10260. 

8271 . 

510. 

711. 

0.0 

5 . 3 

0 . 

1112. 

3414 . 

1470 . 

33.7 

55.9 

9987 . 

2629. 

304 . 

549. 

0.0 

4.0 

0 . 

732. 

3800. 

5516. 

23.4 

36.7 

16239. 

15043 . 

0 . 

829. 

3.6 

4.2 

1657 . 

1396. 

3323. 

5516. 

23.1 

36.7 

16549. 

15243. 

0. 

829. 

3 . 6 

4 . 2 

1716 . 

1396 

2372. 

3099. 

20.4 

28.7 

11627. 

10804 . 

220. 

629. 

0.0 

4 . 3 

979. 

1095. 

4098. 

3933. 

31.7 

27.7 

12927 . 

14379 . 

716. 

961. 

4 . 6 

5 . 9 

1211 . 

1490 . 

4012. 

3451. 

33 

35.4 

12012. 

9754 . 

1069. 

1022. 

0.0 

5 . 4 

1124 . 

1332 . 

4381. 

6233. 

35.0 

35.4 

12517. 

17615. 

1098. 

982. 

0.0 

6 . 3 

1329. 

1835. 

5400. 

6403. 

31.5 

31 . 8 

17142. 

20128. 

0 . 

942. 

6.1 

3.8 

1215 

1502. 

5440 . 

6784. 

29.6 

33.7 

18378. 

20128. 

0 . 

942. 

5.3 

3.8 

1312. 

1599. 

5948. 

6784. 

30 . 1 

33.7 

19760 . 

20128 . 

0 . 

942. 

5.3 

3.8 

1308 

1599. 

5353. 

6784 . 

28.2 

33.7 

18982. 

20120 . 

0 . 

942. 

6.1 

3.8 

1304 

1599. 

4417. 

4332. 

2/ .6 

32.3 

16026. 

13427 . 

706. 

616. 

3.4 

11.0 

2449 

1894 . 

3387. 

4323. 

25.5 

32.3 

13278. 

13399. 

635 . 

614 . 

2.5 

7.8 

1474 

1607 . 

4749. 

4110 . 

30.3 

30 . 9 

15672. 

13297 . 

777 . 

624 . 

1 . 4 

5.7 

2067 

1394 . 

3789. 

3652. 

23.4 

27 . 3 

16216. 

13383 . 

819. 

713. 

3.1 

12.7 

1393 

1976. 

3134. 

3651. 

22.9 

27.3 

10804 . 

13378. 

507. 

713. 

4.4 

10 . 1 

2628 

1740. 

4208. 

3644. 

26.7 

27.3 

14992. 

13352. 

7C5 . 

712. 

3 . 7 

7.8 

1950 

1508. 

4057. 

4497. 

26.2 

32.3 

15382. 

13931. 

836. 

750. 

2.3 

11 . 4 

1953 

1900. 

4138. 

4489. 

30.0 

32.3 

13788. 

13906. 

743. 

749. 

3.0 

8.0 

2373 

1643. 

4746. 

4499. 

29.4 

32.5 

16139. 

13840 . 

777. 

734. 

4 . 5 

6.1 

1621 

1433. 

3825. 

3885. 

31 . 3 

35.0 

12224 . 

11114 . 

640. 

520. 

3 . 1 

12.5 

2818 

1940 . 

2625. 

3883. 

33.8 

35.0 

7765 . 

11108. 

472. 

520. 

2.4 

9.3 

1696 

1672. 

3407. 

3889. 

32.7 

35.3 

10414. 

11027. 

556. 

514. 

3.7 

7.3 

1784 

1475. 

4201 . 

4656. 

36.4 

44 . 7 

11637. 

10412. 

619. 

510 . 

1.6 

10.8 

2264 

1882. 
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Table 7. Continued. 


YIELD 

CRAIN 

UT . 

GRAINS/ SQM 

EARS /SQM 

MAX . 

LAI 

DRY 

MATTER 

0 

P 

0 

p 

0 

p 

0 

p 

0 

p 

0 

P 

3236. 

4452. 

41.6 

42.8 

8250 . 

10408. 

438. 

522. 

1.4 

7.8 

1312. 

1617 . 

3602. 

4388. 

*2.4 

42.8 

8493 . 

10259. 

443 . 

517. 

1.2 

5.2 

0 . 

1366 . 

5529. 

5282. 

32.2 

37.8 

17270. 

13975. 

642. 

510. 

8.8 

11.1 

32*5. 

1931 . 

3793 . 

5274. 

33.0 

37.8 

11493. 

13953. 

536. 

510. 

2.8 

8.2 

1662. 

1663 . 

*397. 

5005. 

33.9 

36.2 

22609. 

13834. 

530. 

518.. 

2.9 

5 7 

1019. 

1*36 . 

*399. 

4506. 

29.7 

30.5 

1*754 . 

15118. 

804. 

5*6. 

4.8 

12.0 

1883. 

1949 . 

3783. 

*370. 

28.9 

29.1 

13090. 

15029. 

605. 

539. 

2.3 

9.0 

1827. 

1691 . 

3655. 

*184. 

31.3 

29.1 

11637. 

14389. 

531. 

532. 

2.0 

6.0 

990 . 

1*32. 

*992. 

57*1 . 

37.4 

41.1 

133*9. 

13969. 

383. 

290. 

0.0 

5.1 

1510 . 

1612 . 

3552. 

57*1 . 

29.3 

41.1 

12124 . 

13969 

348. 

290. 

0.0 

5.1 

1420. 

1556 

6556. 

57*1 . 

43.2 

41 . 1 

15176. 

13969. 

36*. 

290. 

0.0 

5.1 

165 U . 

1612 . 

5436. 

57*1 . 

37.9 

41 . 1 

14344 . 

13969. 

350. 

290 . 

0.0 

5.1 

1300. 

1612 . 

2676. 

4*18. 

22.5 

31 . 6 

11893 . 

13969. 

361 . 

290 . 

0.0 

5.1 

1100. 

1397 . 

4709. 

57*1 . 

36.7 

41 . 1 

12831 . 

13969. 

330. 

290 . 

0 . 0 

5.1 

1690. 

1593 . 

4437. 

5278. 

31.9 

34.7 

13908. 

15229. 

506. 

321. 

0 . 0 

5.4 

11*0 . 

16*5 . 

4386. 

5278. 

27.4 

34.7 

16007. 

15228. 

484 . 

321 . 

0.0 

5.4 

620. 

1594 . 

5427. 

5278. 

32.3 

34.7 

16853 . 

15229. 

511. 

321. 

0.0 

5.4 

6*0 . 

16*5 . 

4862. 

5278 . 

31 . 1 

34.7 

15633 . 

15229. 

4 92. 

321 . 

0 . 0 

5.4 

710 . 

1645 . 

2682. 

3982. 

20.3 

26.1 

13210. 

15228. 

*99. 

321 . 

0.0 

5.4 

620. 

1404 . 

*354 . 

5278. 

30.8 

34.7 

14 X 37. 

15228. 

462. 

321. 

0.0 

* * 

560 . 

1608. 

9108. 

6*17. 

o . u 

41.8 

0 . 

16546. 

0. 

300. 

0.0 

4.2 

0 . 

1562. 

5389. 

6720. 

0.0 

40.6 

0 . 

16566. 

0. 

303 . 

0.0 

4.2 

0 . 

1543 . 

6232. 

7279. 

0.0 

44.0 

0. 

16556. 

0 . 

301. 

0.0 

4.2 

1 . 

1602 . 

8777 . 

7286. 

0.0 

44.0 

0 . 

16573. 

0. 

304 . 

0.0 
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1610 . 

3471 . 

1610. 

0.0 

20.0 

0 . 

9048. 

0. 

265. 

0.0 

3.9 

0. 

9 X 7 . 

4893 . 

4418. 

0.0 

30 . 6 

0 . 

1**17. 

0. 

278. 

0.0 

4.1 

0 . 

1281 . 

6100 : 

5077. 

0.0 

33.5 

0 . 

15176. 

0. 

267. 

0.0 

4 . 0 

0 . 

1348 . 

9158 . 

5092. 

0.0 

33.5 

0 . 

15219. 

0. 

272. 

0.0 

4 . 0 

0 . 

1379. 

2810. 

2682. 

0.0 

31.3 

0 . 

3580. 

0. 

201 . 

0 . 0 

2.9 

0. 

806 . 

5438. 

3641. 

0.0 

33.9 

0. 

10756. 

0. 

207 . 

0.0 

3.0 

0. 

9*1 . 

1389. 

1595 . 

0.0 

20.0 

0. 

7977. 

0. 

202. 

0.0 

2.9 

0. 

681 . 

2281. 

3483. 

0.0 

33 . 9 

0 . 

10290 . 

0. 

183. 

0.0 

2.7 

0 

862. 

2281. 

2043. 

0.0 

30.2 

0 . 

6760. 

0. 

133. 

0.0 

1.7 

0 

550 . 

1*22. 

1713. 

0.0 

25.1 

0. 

6829. 

0. 

139. 

0.0 

1 . 7 

0 

525 . 

1488. 

1786. 

0.0 

28.9 

0 . 

6179. 

0. 

118. 

0.0 

1.4 

0 

481 . 

17. 

316. 

0.0 

26.9 

u . 

Ii 76. 

0. 

85. 

0.0 

0.5 

0 

110 . 

149. 

332. 

0.0 

26.? 

0. 

1237. 

0. 

92. 

0 . 0 

0.5 

0 

116. 

1*9. 

273. 

0.0 

26.9 

0. 

1016. 

0. 

68. 

0.0 

0.4 

0 

93 . 

364. 

281. 

0.0 

26.9 

0. 

1046. 

0. 

71. 

0.0 

0.4 

0 

96. 

4774. 

5818. 

0.0 

45.4 

0. 

12828. 

0. 

194 . 

0.0 

2.5 

0 

1099 . 

4695. 

5978. 

0.0 

45.5 

0. 

131*8. 

0. 

208. 

0 . 0 

2.7 

0 

11*3 . 

4496. 

6162. 

0.0 

*5.6 

0. 

13512. 

0. 

229. 

0.0 

2.9 

0 

1198. 

8133. 

6503. 

0.0 

*9.2 

0 . 

13231. 

0 . 

268. 

0.0 

4.0 

0 

1*67 . 

8926. 

6483. 

0.0 

49.1 

0 . 

13191. 

0. 

263. 

0.0 

3 . 9 

0 

1447 . 

712*. 

6499. 

0.0 

*9.1 

0 . 

13223. 

0 . 

266. 

0.0 

4 . 0 

0 

1460. 

3885. 

6584. 

0.0 

*5. * 

0 . 

14*87. 

0. 

269. 

0.0 

* . 0 

0 

150* . 

6893. 

6635. 

0.0 

45.4 

0. 

14601 . 

0. 

270. 

0.0 

* . 3 

0 

1602. 
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Table 7. Conclnuad. 


YIELD 

CHAIN 

MT. 

CHAINS /SQM 

EARS / 

CQM 

MAX. 

LAI 

DRY 

MATTER 

0 

P 

0 

P 

0 

P 

0 

P 

0 

p 

0 

P 

1736. 

6589. 

0.0 

*5* 

0. 

1**98. 

0 . 

270. 

0.0 

*0 

0 

1511 

2992. 

1223- 

0.0 

20.6 

0. 

156*5. 

0 . 

318. 

0.0 

* . 6 

0 

1307 

63*8. 

6188. 

0.0 

16.6 

0. 

16906. 

0. 

318. 

0.0 

* 6 

0 

162* 

50*8 . 

5***. 

0.0 

12.5 

0. 

16733. 

0 . 

318. 

0.0 

*6 

0 

15** 

*853. 

122*. 

0.0 

20.6 

0. 

156*9. 

0. 

319. 

0.0 

*.6 

0 

1310 

3091. 

1*70. 

0.0 

20.0 

0 

7350. 

0. 

285. 

0.0 

*.8 

0 

9*9 

7108- 

*25*. 

0.0 

37.0 

0 

11*93. 

0. 

283. 

0.0 

* .9 

0 

1276 

355*. 

3391. 

0.0 

29.9 

0. 

113*6. 

0. 

28*. 

0.0 

*8 

0 

11*8 

5091. 

3515. 

0.0 
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Summary 


The enormous diversity of the data base used for testing 
the CERES -Wheat model allowed us to compare every important 
aspect of the model's features. The most important are 
phasic development and developing timing for which almost all 
data sets could be used. The next most important data for 
comparison between the model and experimental results are 
above-ground dry matter production, leaf area development, 
and tillering pattern. We found that measurements of 
tillering were limited, perhaps because tillering occurs 
shortly before or after winter. When more detailed measure- 
ments were taken, they usually were for dry matter partition- 
ing. 


These were the most valuable data sets, because sub- 
routines or parts of subroutines for CERES -Wheat can only be 
validated when their specific output is compared with stem, 
leaf, ear, and root weight measurements. More detailed 
information from experiments would have been valuable. 

Generally, agreement between model output and experi- 
mental results was acceptable or excellent. When 
discrepancies grew too large, usually a conflict could be 
detected, such as heavy diseases or pests in the crops, 
differing definitions of plant parts or units, etc. 

Acceptance of a model's performance remains a personal 
decision and a matter of defining limits. Models restricted 
to certain locations can easily be more accurate when 
regions- ^y applied. CERES, however, was designed for 
universal application- -and it meets the requirements of being 
able to simulate wheat growth and development at any site 
where wheat can be grown. 


CERES -WHEAT NITROGEN 

Some of the important objectives in developing the 
CERES-Wheat-N model were: 

1. To predict response or nonresponse to N 
fertilizer in a diversity of environments. 

2. To predict crop N uptake and N utilization. 

3. To account for the N balance components in 
the soil-crop system. 

4. To predict the time course of biomass 
accumulation and N uptake by the crop. 
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Since model development and testing is somewhat of an 
iterative process in the early stages, most of the data sets 
have been utilized during the development phase and are not 
truly independent. Since development of CERES -Wheat-N has 
necessarily lagged behind the development of CERES-Wheat, the 
opportunity to rigorously test the model with a large base of 
truly independent data sets has not yet arisen. Development 
and testing of the CERES-Maize-N model (Jones and Kiniry 
1986) has proceeded in parallel with the work on the wheat 
model. Because the soil N transformation components of both 
models are identical and since the basic structure of the 
CERES-Uheat model dictates the nature of biomass production, 
yield component determination, and water balance, the testing 
data base can be inferred as having some degree of indepen- 
dence from model development. Thus, in the analyses that 
follow, no attempt has been made to separate truly indepen- 
dent data secs from those used for model development. 

In addition to those test criteria used for testing the 
CERES model, several other procedures were examined: 

1. An approach suggested by Dent and Blackie 
(1979) testing the null hypothesis that the 
intercept coefficients (a) and slope 
coefficient (b) simultaneously are not 
different from zero and unity, respec- 
tively. An F statistic appropriate for 
testing this hypothesis was calculated. 

2. Upper and lower confidence intervals about 
the slope and intercept were also 
determined such that either slope or 
intercept could be identified as signifi- 
cantly departing from the 1:1 line. 

3. A statistic to determine model accuracy as 
defined by Freese (1960) . 

4. A 5% critical error as defined by Reynolds 
(1984). 


Data Base 

Data sets of wheat production from several places in the 
world were assembled for testing and improving the CERES - 
Wheat-N model. Most of the data came from published sources 
and some from unpublished Ph.D. dissertations and other 
unpublished sources. As a complete minimum data set was 
rarely available, the additional climatic and soils 
information was obtained from other reports or personal 
communication. When a few key data were unavailable, certain 


28 



model inputs were estimated using the best available local 
information. 

Comparisons of the predicted and observed time course of 
biomass accumulation of N uptake by the crop and N balance 
components depicted where appropriate observed data were 
available . 

The test data base spans the wheat -growing environments 
from 53 degrees N latitude in the United Kingdom and the 
Netherlands with a 10 -month growing season to the spring 
wheat -growing areas of Canada and the northern United States 
with growing seasons of 90-120 days to the winter-planted 
spring wheat-growing areas of the Middle East and Australia. 
A diversity of soil types and fertilizer application 
patterns, sources, and timings is also represented in the 
dati base. 


Individual Data Sets 

Garden City. Kansas. U.S.A. (1980 and 1981) . The 
experimental design was six N rates (0, 28, 56, 84, 112, and 
140 kg N/ha) with four irrigation timing strategies (preplant 
irrigation only, preplant + irrigation at jointing, preplant 
+ irrigation at flowering, irrigation at all three times). 
The soil was a clay loam (Aridic Argiustoll) and fertilizers 
were applied broadcast at preplanting followed by incorpor- 
ation. The variety was Newton. The experiment was conducted 
by Dr. Mark Hooker; data were obtained by personal communi- 
cation and are reported in Wagger (1983). 

Manhattan. Kansas. U.S.A. (1981) . The experimental 
design was three N rates (0, 60, and 180 kg N/ha) and a plus 
and minus irrigation treatment. The 180 kg N/ha rate was 
divided; half was applied at planting and the other half at x 
days after planting (Wagger 1983). Additional data were 
obtained by personal communication (Drs. M. Wagger and D. 
Kissel, Kansas State University). 

Hutchinson. Kansas. U.S.A. (1979 and 1980) . The 
experimental design was six N rates (0, 28, 56, 84, 112, and 
140 kg N/ha) applied preplant followed by incorporation. In 
1979, the variety was Centurk; in 1980, Newton was used. 

Swift Current. Saskatchewan. Canada (1975) . The 
experiment comprised seven N rates (0, 20.5, 41, 61.5, 82, 
123, and 164 kg N/ha) with a plus and a minus irrigation 
treatment (Campbell et al. 1977 a,b). Following planting, 15 
cm diameter lysimeters were driven into the soil to a depth 
of 120 cm. Five harvests during the growing season were 
made. The variety was Manitou. Climate and soils data were 
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obtained from the authors. The instrument gathering solar 
radiation data malfunctioned for three weeks during the early 
grain filling stage. The missing data were estimated by 
fitting a function to radiation of maximum and minimum 
temperature, presence or absence of rainfall, and the day of 
the year. The reliability of these estimates is not known. 

Northwest Syria (1979 and 1980) . Experiments were 
conducted at four sites in Aleppo province by Dr. M. Stapper. 
At three of the sites (Brida, Jindiress, and Kafr Antoon) , 
two N rates (0 and 60 kg N/ha) were applied. At each site, 
three spring wheat varieties (Hexipak, Sonalika, and Novi 
Sad) were compared. At the fourth site, additional 
irrigation treatments were added. At Kafr Antoon, a late 
frost was suspected, and at Jindiress, the variety Mexipak 
suffered from rust (Dr. H. C. Harris, personal communication, 
ICARDA) . Data are reported in Stapper (1984) . 

Wonean Hills. Western Australia (1966) . The experiment, 
using the spring wheat variety Gamenya, was designed to 
examine the fate of anhydrous ammonia and urea applied to a 
loamy sand (Mason and Rowley 1969) . Since the model 
simulates both of these fertilizer materials as identical 
ammoniacal sources, comparisons were made with the mean of 
these two treatments. The experiment showed no significant 
differences between the sources. The rates of N applied were 
0 and 61 kg N/ha applied preplant. Solar radiation data were 
estimated from recorded hours of sunshine. 

Lancelin. Western Australia (1967) . The experiment was 
on a very coarse, siliceous sand and was designed to examine 
the fate of urea applied at various intervals after planting 
(Mason et al. 1972). Urea was applied at 77 kg N/ha either 
at planting or at 2, 4, or 8 weeks after planting. Delaying 
the application resulted in an almost threefold increase in 
grain yield. Solar radiation for Ferth, 100 km distant, was 
used as part of the climatic data. 

Rothamsted. England (1975) . The experiment compared 
three varieties of winter wheat (Maris Huntsman, Capelle 
Desprez, and Maris Fundin) over either rates of N (0, 30, 60, 
90, 120, 150, 180, and 210 kg N/ha) (Pearman et al. 1978). 
The fertilizer was applied 163 days after planting. Few 
significant variety X N interactions were recorded. The 
initial soil mineral N values and soil water contents were 
interpolated from estimates provided by the authors. 

Flevopolder. The Netherlands (1975) . The experiment was 
designed to test the effects of late applications of N on 
leaf area duration, assimilation nutrient uptake, and growth 
of grains (Spiertz and Ellen 1978) . A zero N treatment was 
not included, thus comparisons of apparent fertilizer 
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recovery were not made. The soil was a marine clay in a 
reclaimed polder; the winter wheat variety was Lely. A water 
table was present during the course of the growing season. 
This was simulated by assuming that layers in the profile 
below 1 m deep were filled to saturation. Without this 
assumption, the model predicted considerable moisture stress. 
Soil mineral N analyses were made after the crop was planted, 
thus estimates based on these were used for the initial 
mineral N input values supplied to the model. 

Wageningen. The Netherlands (1977) . These data are 
reported by Ellen and Spiertz (1980). The experiment 
examined various strategies of splitting fertilizer 
applications on uptake and yield of grain. The variety used 
was the winter wheat Donata. The soil was repotted as a 
river clay with 45% silt. 

Carrington. North Dakota. U.S.A. (1969-73) . These data 
were reported by Bauer (1980). The experiments involved a 
comparison of several varieties of hard, red spring wheats 
over several rates of nitrogen and over five years. Half of 
the experiment was irrigated and the remainder dryland. 
Since this experiment yielded a massive data set, only 
selected contrasting years were utilized in the testing data 
base to avoid biasing the data base with to many points from 
one location. Straw yields and .straw N percent were not 
reported in some instances and were estimated from the grain 
yields and grain protein concentrations using the regression 
procedures described in the publication. This may lead to 
some errors in estimation of observed biomass and N uptake . 

Wagga Wagga. N.S.W.. Australia (1962) . The experiment, 
reported by Storrier (1966), used four rates of N. Half of 
the experiment was irrigated; the other half was dryland. As 
a split plot treatment, a later application of 45 kg N/ha as 
sodium nitrate was made to half of the plots. Storrier 
(1966) reports a negative response to applied N, but 
examination of the variation in the experiment indicates more 
of a case of nonresponse to N. Some lodging was reported in 
the high N plots, which the model would not have been able to 
account for. The initial mineral N in this experiment was 
very high. Several gaps existed in the weather record for 
this experiment. Solar radiation was estimated from hours of 
sunshine data recorded from the site or nearby. Where the 
record was blank for both sites, mean values for that time of 
year were used. 
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Results and Discussion 


Validation: Difference Measures and Summary Statistics 

Grain Yield . Simulated grain yields are tabulated 
against observed counterparts (Table 10) . The means and 
standard error of predictions closely approached those of the 
observations (Table 8). The degree of scatter around the 1:1 
line (Figure 8a) is very small. Thirty yield predictions 
from 240 deviated more than one standard deviation from their 
observed counterparts. Data sets, where predictions were 
poor, were those from Carrington, North Dakota, U.S.A. 
(1969), Hutchinson, Kansas, U.S.A. (1980), Flevopolder, The 
Netherlands (1975), Jindiress, Syria (1980), and some 
individual treatments from some of the remaining data sets. 

The regression line (Figure 8a) has a slope greater than 
unity (1.033) which differs from unity by slightly more than 
5% confidence interval (0.030). The F statistic for the 
regression is significant at the 5% level due to the slope. 
This slope of greater than unity and the small positive 
intercept (145.6 kg/ha) implies that the model has some 
tendency to overpredict yields across the range of 
observations. This is further indicated by the small 
positive MBE (Table 8) . This finding is not surprising and 
it is not sufficient cause to reject the model since the 
assumption is made that all nutrients other than N were 
present in nonlimiting quantities and that other factors not 
accounted for by the model (pests and diseases, crop lodging 
effects, frost induced sterility, etc.) had no influence on 
yield. As mentioned in the description of the data sets, 
these assumptions may not always have been entirely 
fulfilled. No quantitative data were available to suggest 
that some of the test data sets should have been eliminated 
from the testing data base . 

All other statistical criteria for model evaluation 
(correlation coefficient, chi-square test, and modified 
Freese statistic) indicated the model was acceptable. The 
critical error term of Reynolds (1984) indicates that the 
model will predict grain yield within an error of 1,865 kg. ha 
with a 95% confidence. 

Biomass . Simulated biomass is tabulated against its 
observed counterparts (Table 10) . Means and standard error 
for predicted biomass closely resemble those from the 
observed data (Table 8). More scatter about the 1:1 line 
occurred for biomass predictions (Figure 8b) than for grain 
yield. The greater spread is generally from a poorer 
simulation of straw yield. Most of the 55 points falling 
outside the bounds of the + 1.0 standard deviation were from 
whole data sets rather than from individual treatments across 


32 



a range of data sets. Data sets where simulation of biomass 
was generally poor were Wageningen, 1977; the variety 
Nugaines at Pendleton and Dufur, Oregon, 1971; various 
treatments within the Jindiress, 1980; Tel Hadya, 1980; Kafr 
Antoon, 1980; and Garden City, 1981. The simulated biomass 
for the Hutchinson 1980 data set was consistently low. 

The slope of the regression line significantly deviated 
from the 1:1 line indicating a tendency of the model to 
overpredict biomass particularly at the high end of the 
range. This is also indicated by an MBE of 1,406 kg/ha. The 
correlation coefficient, chi-square, and modified Freese 
statistic all indicated model predictions were acceptable. 

Total N Uptake. Grain N Uptake, and Grain Protein . 
Performance of the model in predicting these parameters was 
generally poorer than the simulation of grain yield. Forty- 
four points from a total of 223 fell outside of the bounds of 
± standard deviation of the 1:1 line for total N uptake for 
grain protein and grain N uptake: 31 from 215 and 60 from 
215, respectively. The slope of the regression line (1.042) 
for total N uptake was just beyond the 5% interval for slope 
(1.039). Similarly, the intercept (7.0) lies beyond the 5% 
confidence interval for the intercept (6.34). There is a 
tendency for the simulations to exceed the observations 
although the correlation coefficient, chi-square test, • and 
modified Freese procedure all indicate the simulations are 
acceptable . 

Total N uptake was consistently underpredicted for the 
Swift Current, Canada, data, and consistently overpredicted 
for the Waite Institute, Dufur, and Pendleton data sets. 
Grain N uptake was simulated fairly closely (Figure B.l). 
Much of the error involved in simulation of total N uptake 
arose from poor simulation of the concentration of N in the 
straw at harvest. 

The range of simulation values was consistently less 
than that observed (0.1 to 1.3%). Some of these differences 
may occur due to differences in harvesting technique and time 
of harvest. If significant amounts of chaff or leaf 
materials are not included in the sample, the reported straw 
N concentration also will be low. The model makes no attempt 
to account for losses of N from the vegetative material 
through leaching of N compounds from harvest ripe straw or 
via volatile losses from senescing leaves. Several data sets 
had less N in plant top tissue at harvest than at anthesis, 
indicating some losses . 

The scatter of points around the 1:1 line (Figure 8) was 
much higher for grain protein (grain N percent multiplied by 
5.7) than for many other parameters. Both slope and 
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intercept of the regression line are significantly different 
from the 1:1 line. The chi-square test also indicates the 
simulations are significantly different from the observ- 
ations, but the modified Freese statistic indicates the model 
is still acceptable. 

The simulation of grain protein concentration has been 
to date one of the most difficult components in the whole 
model to get working satisfactorily. In several of the data 
sets, grain protein concentration was consistently overpre- 
dicted or underpredicted. It was difficult to determine if 
in any of these cases a genotypic factor was involved. 
Adding a further genetic coefficient to the model input data 
requirements to help explain cultivaral differences in grt-in 
protein accumulation has so far been avoided. Further 
investigation of this aspect of the model is warranted. 

Drv Weight and N Uptake an Anthesis . In many of the 
studies, harvests were made at or near anthesis. The 
simulated data used for the comparisons were the correspond- 
ing values for N uptake and biomass on the date of harvest 
(i.e., not necessarily on the simulated date of anthesis). 
Biomass was generally overestimated at anthesis. The 
regression line significantly deviates from the 1:1 line 
(Figure 8g) and 72 points of the 161 fell outside the bounds 
of + 1.0 standard deviation of the 1:1 line. The chi-square 
test indicated the predicted biomass differed significantly 
from that observed. The simulated N uptake at anthesis 
showed much less scatter than the predictions for biomass. 
The simulations were acceptable within all of the statistical 
criteria examined. The model substantially underestimated 
the anthesis N uptake of several of the treatments from the 
experiments at Wagga Wagga, 1962. The resulting slope of the 
regression line thus is a little less than 1.0. Some 
compensating errors in the simulation of plant N concen- 
tration may occur if the simulated biomass is incorrect and 
the simulated N uptake is correct. The MBE terms (Table 8) 
also indicate a large overestimate of biomass at anthesis and 
a slight underestimate of N uptake at anthesis. 

Kernel Weight and Kernels Per Square Meter . Overall, 
the model had a slight tendency to underestimate the number 
of grains per square meter and overestimate the weight of 
individual kernels. The model consistently overestimated 
kernel weight for the variety Capelle Desprez in the 
Rothamsted 1975 experiments. For the final determination of 
grain yield, small errors in either of these components ire 
of no consequence provided compensation occurs (i.e. , low 
kernel number is compensated by a high kernel weight) . These 
two yield components are, however, important indications to 
timing of certain stresses, and for the model to be useful, 
they should be reasonably correct. For both parameters, the 
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slope of the regression line significantly differs from the 
1:1 line, but the intercepts were within the confidence 
interval. There was a noticeable tendency for the model not 
to display the same sensitivity in kernel weight to rates of 
applied N as the observed data did. 


Apparent Recovery . Apparent recovery (AR) is a 
parameter often used in fertilizer research to indicate the 
efficiency of fertilizer use. From a modelling standpoint, 
it is a particularly challenging parameter on which to test 
the model, since AR depends on the accurate simulation of two 
treatments simultaneously. It is calculated: 


AR - 


NliP _ - NUP 
f u 

Rate 


x 100 


where NUPf - N uptake from a fertilized treatment 

NUP U - N uptake from an unfertilized treatment 
Rate - rate of fertilizer applied. 

Since AR depends on the N uptake from two different 
treatments, small errors in the prediction of either can lead 
to quite spurious values for the calculated AR. This is 
indicated in Table 9 where 10% errors in prediction of both 
the fertilized and unfertilized treatment leads to errors in 
the calculated AR of 33% and 38%. 


Table 9. Effect of Errors in N Uptake Simulation on Errors and 
Apparent Recovery. 



NUP f 

NUP U 

Race 

AR 

Error (Z) 

Sample Obs 

Z0 

40 

30 

67 

- 

10Z Error 1 

66 

36 

30 

100 

33 

10Z Error 2 

5* 

44 

30 

33 

-3* 


The means and standard errors of the predictions closely 
approximated those of the observed. Most of the statistical 
parameters (Table 8) indicate a significant difference 
between predicted and observed values . The modif iees to N in 
both years at Garden City, Kansas, were not apparent in any 
of the irrigation treatments. The simulations reasonably 
approximated the observations across the range of treatments. 

Predicted N responses for the Syrian data sets were 
generally very good. Yields for some of the varieties were 
overestimated at Jindiress and Kafr Antoon, but the 
simulations were consistent with the remarks noted in the 
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above section "Description and Testing Data Base." Yields 
for the longer duration variety Novi Sad were overestimated 
at Tel Hadya, 1979, and small underestimates of yield 
occurred for the zero N treatments at Brida . 

Responses to N applied either at planting or in split 
applications ;ere not apparent at Madras, Oregon, or for 
single applications at Dufur. Simulations for these treat- 
ments were very close to the observations. The model 
underestimated yields for the variety Hyslop at Pendleton, 
but overestimated yields across the range of N rates for the 
variety Nugaines . 

The model underestimated the grain yield at low N rates 
in the Swift Current, 1975, experiment, but simulations for 
the remainder of the response curve were excellent. 
Excellent yield simulations were also recorded for the 
Rothamsted data sets. 

The sensitivity of the model to differences in 
fertilizer application pattern (timing) is well illustrated 
by the data sets from Lancelin, Wageningen, and Bozeman. 

Grain protein concentration generally was well simulated 
across the range of N rates, except for the cases noted above 
(section "Total N Uptake, Grain N Uptake, and Grain 
Protein"). These exceptions rendered the simulation of grain 
protein to be unacceptable statistically when all data sets 
were combined. Given the constraints noted above in model 
development and the proximity of most of the simulations 
depicted in Figure B.6, there is insufficient evidence to 
reject the model for the applications for which it was 
designed. 

The model tended to overestimate total N uptake, but 
there was no consistent pattern of overestimation or under- 
estimation. Despite these problems, the model captures most 
of the observed effects of N uptake for most of the data 
sets. Some further study is required to elucidate the 
problems with the North Dakota and Syrian data sets. 


Seasonal Patterns of Biomass and N Uptake Accumulation 

Across the range of data sets, predicted biomass was 
slightly out of phase with observed biomass. The model 
tended to predict higher biomass accumulations earlier in the 
season than the observations would indicate. While there was 
a noticeable trend, the errors were not large and were 
consistent with those observed for the non-nitrogen version. 
Errors were large, however, at the low N rates for the Swift 
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Curient data set and for some treatments in the Tel Hadya 
1980 data set. 

A similar pattern in seasonal N uptake to that of 
biomass was observed (early overestimation of N uptake) . 
Seasonal patterns of N uptake were poorly simulated at the 
low N rates in the Swift Current data set, but reasonably 
simulated in the Kansas data sets. 


Seasonal Patterns of N Balance 

When attempting the validation of the soil N components 
of the model, it was originally intended to attempt a layer- 
by-layer comparison of each of the predicted nitrate and 
ammonium concentrations with those observed. Analysis of the 
observed data in most instances first indicated very large 
standard errors and other seeming anomalies were sometimes 
apparent. In some of the data sets, mineral N concentrations 
were low after fertilizer addition and increased as the crop 
grew and withdrew N from the soil. While this may be 
indicative of turnover occurring within the soil systems due 
to microbial activity or ammonium adsorption/desorption on 
clay surfaces, the anomalies in the individual layers were so 
gross as to discard many of them. The errors associated with 
the layer data for most data sets were such that the 
simulations may have been in error by some 200% and yet still 
be within the error bounds of the observatir s. 

To provide some meaningful validation of these 
components of the model, total mineral N in the soil at 
various times was used as a test criterion. While many of 
the error conditions noted above will still affect this total 
soil pool, the values obtained were more consistent with what 
could be reasonably expected. 

Across the range of data sets studied, some within- 
season differences in predicted and observed balances occur, 
but given the magnitude of the errors cited, the N balance 
simulations are plausible. 
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Figure 8. Comparison of Predictions of the CERES -Wheat-N Model with 
Observed Data from Field Experiments. 

a. GRAIN YIELD b. BIOMASS 









Table 8. Stannary Measures -for ell Data Sats. 


Varieble 

Units 

N 

6 

P 

SO 

SP 

a 

b 

R 

D 

F 

Biomass 

kg/ha 

222 

10.313 

11.719 

3.375 

3.897 

189.3 

1.118 

0.82 

0.86 

15.907 

Grain yield 

kg/ha 

210 

3,953 

*.227 

1.716 

1.719 

1*5.6 

1.033 

0.8* 

0.9) 

3.263 

Total N uptake 

kg N/ha 

223 

110 

121 

*7 

*5 

7.0 

1.0*2 

0.72 

0.83 

7.8*38 

Gram N uptake 

kg N/ha 

21 5 

84 

87 

3* 

31 

4.4 

0.996 

0.7* 

0.85 

1.298 

Grain protein 

Z 

215 

12.3 

’1.7 

2.95 

3.55 

0.198 

0.939 

0.55 

0.74 

2.933 

An thesis DW 

kg/ha 

161 

7.25* 

8.681 

2.166 

2,818 

261.7 

1.161 

0.58 

0.70 

22.771 

An thesis A uptake kg N/ha 

151 

118 

116 

57.1 

39.31 

11.66 

0.887 

0.67 

0.78 

2.3370 

GPbii 

NO 

152 

12.381 

11.861 

*.668 

5.32* 

362 

0.929 

0.63 

0.78 

1.080 

Kamel weight 

mg 

13* 

33.6 

37.9 

7.02 

3 . *0 

0.99E 

1.096 

0.3* 

0.59 

18.69 

App recovery 

Z 

137 

*5.7 

*3.* 

29* 

25.68 

8.37 

0. 7 

0.37 

0.65 

6.2* 


Variable 

Units 

N 

P 

CIS’ OP 

CIINT 

* 2 

MAE 

. MBE 

RMSE 

RT 

£ o 

Biomass 

kg/ha 

222 

0 

0.028 

*23.6 

136.01 

2. 2*8.1 

1. *05. 590 

2.6*8.19 

13.06 

*ao 9 

Grain yield 

kg/ha 

2*0 

0.0*00 

0.030 

173.9 

85.08 

806.08 

27*. *3 

1.02*. 16 

15.39 

156.5 

Total N uptake 

kg N/ha 

223 

0.0005 

0.0392 

6.3* 

133.03 

28.55 

11.535 

36.27 

22.43 

65.9 

Gram N uptake 

kg N/ha 

215 

0.2752 

0.035* 

4.344 

105.13 

19.* 

3.491 

2*. 0*0 

13.64 

43.6 

Gram protein 

Z 

215 

0.055* 

0.0330 

0.582 

250.832 

2.363 

-0.550 

3.192 

18.65 

5.3 

An thesis DW 

kg/ha 

161 

0 

0.0509 

533.22 

259.72 

2.159.6 

1.426.76 

2.759.73 

19.26 

*9* 6 

Anthesis N uptake 

kg N/ha 

151 

0 . 1002 

0.0*67 

8.239 

82. *6 

3*. 36 

-1.660 

*2.352 

23.62 

75 7 

GPSM 

NO 

152 

C.3*22 

0.0516 

362.30 

131. *3 

3. *62. 3 

-520.78 

*.355.36 

22. *3 

778.6 

Kernel weight 

mg 

1*0 

0 

0.0*52 

2.172* 

27*. 77 

8.091 

4.2286 

9.3738 

12.00 

17.5 

App recovery 

Z 

137 

0.0026 

0.0850 

6.032 

152.20 

23. *60 

-2.281 

31. 177 

69. *8 

55.5 


N * Humber of observations. 

0 ■ Man of observations . 

P * Mean of predictions. 

SO ■ Standard deviation of observations. 

SP « Standard deviation of predictions. 

a ■ Intercept term from regression of predicted on observed, 
b “ Slope term from regression of predicted on observed. 

R ■ Regression coefficient. 

D « Index of agreement (Willmott 1982). 

F - F statistic calculated as per Dent and Blackie (1979). 

P * Probability of exceeding F. 

CTSLOP * SZ confidence interval about slope of regression line. 
CIINT * SZ confidence interval about intercept of regression line. 
X i “ Chi-squar* 

MAE * Mean absolute eiror (Willmont 1982). 

MBE - Mean bias error ((•illmont 1982). 

RMSE • Root mean square error. 

RT " Model accuracy (Freese 1960). 

E* ■ SZ critical error as defined by Reynolds (1981). 
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Tabl« 10. Listing of N-Modal Tasting Data Basa (0 « Obsarvad Valua and P “ Pradletad Value). 


CR/ IN 

YIELO 

BI0KAS5 

N UPTAKE 

CRAIN N 

NUPTAKE 

GRAIN 

PROTEIN 

ANTHESIS 

N UPTAKE 

0 

P 

0 

P 

0 

P 

0 

P 

O 

P 

0 

P 

406*. 

1727. 

9210. 

9518. 

72.9 

102.7 

611 

78.5 

8.6 

12.0 

58.6 

106 1 

5129. 

1921. 

11*75. 

10126 . 

122.1 

121.0 

98.2 

82.6 

10.5 

12.0 

97.3 

i:-6.7 

5126. 

*201. 

1101*. 

11015. 

127.8 

1*2.9 

104.1 

88* 

11.6 

12.0 

92.5 

1*7.0 

19*9. 

5090. 

8911. 

121*9. 

61.1 

115.1 

51.1 

107.2 

7.* 

12.0 

62.* 

138 9 

5080. 

5521. 

11772. 

117*0. 

10*. 1 

16*. 2 

81.8 

116.2 

9.* 

12.0 

84.8 

169.2 

5082 

6011. 

11916. 

1518* . 

127.5 

218.8 

101.4 

126.5 

11.6 

12.0 

116.3 

22* . 9 

*62*. 

*118. 

10118. 

10*02. 

91.2 

101.6 

75.* 

81.1 

9.3 

11.5 

72.8 

106.1 

5018. 

*160. 

1108*. 

11295. 

116.5 

122.5 

91.6 

91.8 

10.4 

12.0 

87.3 

126 7 

*9*8. 

*575. 

11127. 

11981. 

119.6 

1*1.0 

96.* 

96.1 

11.1 

12.0 

15.3 

1*7 0 

*175. 

5090. 

9186. 

121*9. 

76.1 

11*. 7 

61.0 

107.2 

8.6 

12.0 

76.6 

138.9 

5211- 

5521. 

12187 . 

117*0. 

115.0 

16*. 2 

92.1 

116.2 

10.1 

12.0 

98.1 

169.2 

*988. 

6011. 

11916. 

15186 . 

129.9 

218.8 

101.5 

126.5 

11.6 

12.0 

109.9 

22* 9 

*585. 

1799. 

10005. 

9776. 

90.7 

106.5 

74.0 

80.0 

9.2 

12.0 

0.0 

0.0 

5185. 

1915. 

11216. 

1028* . 

121.8 

121.6 

99.2 

82.8 

10.9 

12.0 

0.0 

0.0 

5*00. 

*010. 

12*16. 

10521. 

118.9 

129.2 

111.8 

84.8 

11.8 

12.0 

0.0 

0.0 

5119. 

5269. 

11501. 

12811 . 

97.7 

1*4.6 

80.8 

110.9 

9.0 

12.0 

0.0 

0.0 

5185. 

576*. 

12897. 

1*287. 

125.7 

181.9 

95.5 

121.1 

10.5 

12.0 

0.0 

0.0 

5**0. 

5923- 

11118. 

1*851. 

115.9 

202.5 

110.7 

124.7 

11.6 

12.0 

0.0 

0.0 

*756. 

*208. 

10*27. 

1061*. 

101.* 

106.2 

8*. 1 

8*. 9 

10.1 

11.5 

0.0 

0.0 

5188. 

*181. 

11*12. 

11203. 

12*. 5 

121.0 

100.1 

92.2 

11.0 

12.0 

0.0 

0.0 

59*8. 

**57. 

11079. 

11*18. 

1*8.* 

129.2 

117.9 

91.8 

11.1 

12.0 

0.0 

0.0 

5079. 

5269. 

1190*. 

12811 . 

100.0 

1**.6 

82.9 

110.9 

9-1 

12.0 

0.0 

0.0 

5109. 

57**. 

12187. 

1*287. 

122.6 

181-9 

96.8 

121.1 

10.8 

12.0 

0.0 

0.0 

5151. 

5921. 

110*1. 

1*851 . 

117.7 

202.5 

108.9 

12*. 7 

11.6 

12.0 

0.0 

0.0 

2117. 

2226. 

599*. 

*96*. 

59.2 

*6.1 

**.l 

17.9 

10.9 

9.7 

50.9 

* 6.5 

1110. 

1011. 

10178. 

72*5. 

101.2 

72.8 

61.7 

58.6 

10.9 

11.1 

101.0 

7*1 

*521. 

1717. 

126*9. 

9088. 

166.6 

112.8 

99.9 

7*.l 

12.6 

11.* 

175.0 

115.7 

1*18- 

16*6. 

1926. 

1871. 

17.7 

11.1 

26.7 

2*. 5 

10.6 

8.5 

*5.5 

32.9 

1025. 

2769. 

9*2*. 

6760. 

87.8 

61.1 

56.8 

*9.1 

10.7 

10.1 

81.5 

66.3 

*695. 

*168. 

1106*. 

10689. 

166.5 

111.0 

95.5 

81.* 

11.6 

11* 

1*7.6 

133.9 

1828- 

2075. 

61*7. 

5215. 

66.6 

58.7 

18.5 

50.2 

12.0 

13.8 

58-9 

61.7 

222*. 

2179. 

8918. 

6165. 

97.1 

72.7 

50.7 

62.2 

11.0 

1*.9 

86.7 

76.2 

?*!75. 

2577. 

7*81. 

7095. 

79.6 

89.* 

52.2 

76.9 

12.8 

17.0 

8*. 1 

93.7 

216*. 

2759. 

8175. 

8101. 

98.1 

111.6 

50.5 

89.5 

13.3 

18.5 

112.8 

117.6 

2165. 

2958. 

9202. 

9026. 

107.5 

129.1 

59.1 

99.6 

1*3 

19.2 

122 . 2 

136.2 

2172. 

1079. 

95*0. 

9**0. 

119.7 

119.0 

58.7 

101.7 

l*.l 

19.2 

100.4 

1*6.6 

Ijj.7 . 

158* . 

9176. 

5511. 

80.1 

70.8 

67.1 

17.2 

10.9 

13.4 

0.0 

73.1 

*060. 

2021 . 

11225. 

6*90. 

107.7 

91.1 

86.9 

*7.5 

12.2 

13.* 

115.1 

84.0 

1196. 

211*. 

886*. 

7086. 

92.* 

90.* 

79.2 

5*. 9 

13.3 

13.4 

117.1 

92.5 

1999. 

2*95. 

11221. 

7*22. 

119.1 

96.1 

9*. 7 

58.7 

13.5 

13.* 

13*. 0 

98.6 

1786. 

2579. 

10289. 

7598. 

117.1 

95.9 

91.7 

60.6 

l*.l 

13.4 

133.8 

102.3 

1*17. 

2628. 

9205. 

7699. 

105.* 

101.7 

86.9 

61.8 

1*.5 

13.4 

114.8 

10*. 2 

1819. 

2*22. 

10592. 

6987. 

81.6 

86.1 

70.0 

56.9 

10.4 

13.* 

130.7 

88* 

*1*0. 

1201. 

11*51. 

86*1 

95.1 

101.5 

78.* 

75.1 

10.8 

13.* 

121.9 

106.1 

**15. 

1812. 

12727. 

100C 6. 

101.2 

122.0 

89.1 

89.6 

11.5 

13.* 

1*7.9 

12*. 8 

4168. 

*198. 

11278. 

10911 . 

121.0 

114.1 

9*. 1 

98.7 

12.3 

13.* 

152.2 

137.1 

'•?BB 

*155. 

11C34. 

11298. 

1*7.0 

119.1 

101.1 

102.* 

13.7 

13.* 

176.6 

1*2.7 

1819. 

***S. 

1169*. 

11525. 

11*. 9 

1*2.1 

90.* 

104.5 

11.5 

13.* 

154.7 

1*5.5 
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Table 10. Continued. 


CRAIN 

YIELD 

BIOMASS 

N UPTAKE 

GRAIN N 

HUPTAXE 

GRAIN 

PROTEIN 

ANTHESIS 

N UPTAKE 

0 

P 

0 

P 

O 

p 

0 

p 

0 

P 

0 

p 

3464. 

1593. 

9324. 

5531. 

80.9 

70.8 

67.5 

37.4 

n.i 

13.4 

132.2 

73.3 

3789. 

2021. 

10079. 

6430. 

101.3 

81.4 

82 4 

47.5 

12.4 

13.4 

115.1 

84.0 

3685 . 

2336. 

9734. 

7111. 

99.7 

90.4 

83.4 

54.9 

12.9 

13.4 

117.1 

92.5 

3819. 

2497. 

10815. 

7480. 

116.3 

96.3 

91.8 

58.7 

13.7 

13.4 

134.0 

98.6 

3886. 

2583. 

10501. 

7685. 

127.3 

99.9 

98.2 

60.7 

14.4 

13.4 

133.8 

102.5 

3678. 

2629. 

10148. 

7798. 

113.6 

101.8 

91.0 

61.8 

14.1 

13.4 

114.8 

104.2 

3946. 

2424. 

11047. 

6991. 

100.2 

86.1 

80.3 

57.0 

11.6 

13.4 

130.7 

88.4 

4221. 

3201. 

12837. 

8644. 

115.0 

103.5 

87.4 

75.3 

11.8 

13.4 

123.9 

106.1 

4288. 

3812. 

13448. 

10006- 

117.3 

122.0 

92.5 

89 6 

12.3 

13.4 

147 9 

124.8 

3959. 

4198. 

4791. 

10913. 

132.1 

134.3 

91.0 

98.7 

13.1 

13.4 

152.2 

137.1 

4026. 

4355. 

1. *0. 

11298. 

134.5 

139.1 

99.6 

102.4 

14.1 

13.4 

176.6 

142 . 7 

3933. 

4445. 

13723. 

11525. 

131.4 

142.1 

95.2 

104.5 

13.8 

13.4 

154.7 

1*5. 5 

1617. 

788. 

3781. 

1398. 

52.9 

15.9 

44.0 

14.4 

15.5 

10.4 

35.8 

11.8 

1578. 

1279. 

3957. 

2754. 

53.8 

24.9 

43.5 

20.9 

15.7 

9.3 

42.1 

2*.l 

1754. 

1528. 

4603. 

3903. 

65.7 

36.2 

51.1 

29.5 

16.6 

11.0 

49.6 

36.6 

1732. 

1832. 

5042. 

5016. 

70.2 

45.8 

51.4 

36.6 

16.9 

11.4 

59.0 

51.7 

2028. 

2067. 

5063. 

5842. 

78.9 

63.3 

61.9 

41.7 

17.4 

11.5 

58.9 

65.2 

2192. 

2332. 

5151. 

6497. 

81.9 

74.5 

64.6 

47.0 

16.8 

11.5 

65-6 

77.1 

2367. 

2465. 

J'i40. 

6776. 

68.7 

80.9 

70.6 

49.7 

17.0 

11.5 

60.1 

83.6 

2754. 

602. 

6121. 

1096. 

82.3 

11.6 

68.1 

10.4 

14.1 

9.8 

46.4 

9.5 

3792. 

1872. 

8494. 

3406. 

116.7 

36.0 

96.5 

32.2 

14.5 

9.8 

50.6 

32.2 

4154. 

2875. 

9031- 

5297. 

134.6 

55.5 

110.0 

49.4 

15.1 

9.8 

70.5 

50 * 

4395. 

3906. 

9431- 

7295. 

142.5 

74.9 

114.1 

66.5 

14.8 

9.7 

94.7 

70 6 

4395. 

4583. 

10122. 

8911. 

148.5 

91.0 

116.4 

80.4 

15.1 

10.0 

105.3 

90.0 

5425. 

5305. 

994o. 

10435. 

170.5 

111.3 

14* . 7 

97.7 

15.2 

10.5 

108-2 

113.3 

4883. 

5523. 

10142. 

10794 . 

165.8 

118.8 

134.5 

104.6 

15.7 

10.8 

116.8 

12*2 

1420. 

949. 

4735. 

3853. 

48.9 

52.0 

32.1 

21.3 

12.9 

12.8 

46.5 

54.0 

1430. 

1327. 

5010. 

5768. 

62.3 

97.1 

40.4 

29.8 

16.1 

12.8 

48.9 

100.5 

1940. 

1223. 

5430. 

4463. 

57.6 

51.6 

42.5 

32.4 

12.5 

15.1 

38.0 

53.3 

1910. 

1811. 

S430. 

6723. 

64.2 

96.0 

47.2 

48.0 

14.1 

15.1 

50.2 

99.9 

1360. 

1100. 

5200. 

3898. 

44.0 

51.3 

26.5 

39.8 

11.1 

20.6 

39.9 

54.3 

1340. 

1434. 

5107. 

5730. 

61.4 

97.1 

36.7 

51.8 

15.6 

20.6 

51.2 

101.0 

2510- 

4735. 

7000. 

11537. 

49.4 

113.9 

37.9 

92.2 

8.6 

11.1 

39.8 

111.6 

3950- 

5511. 

10925. 

14101. 

106.9 

155.8 

75.5 

107.3 

10.9 

11.1 

110.7 

156.5 

2840. 

5291. 

7350. 

10495. 

46.4 

104.8 

38.9 

91.0 

7.8 

9.8 

62.5 

*9.8 

4330. 

6270. 

11160. 

13806. 

101.1 

151.6 

76.7 

132.0 

10.1 

12.0 

120.2 

95.5 

■•070. 

3617. 

5816. 

10970 . 

40.2 

115.0 

28.0 

95.2 

7.7 

15.0 

34.2 

117.7 

3890. 

4258. 

10167. 

13586. 

82.7 

164.0 

63.5 

120.3 

9.3 

16.1 

86.5 

167.5 

4410. 

3151. 

12710. 

8241. 

0.0 

88.8 

0.0 

75.7 

0.0 

13.7 

0.0 

91.1 

4090. 

4112. 

12560. 

11246. 

0.0 

133.1 

0.0 

102.4 

0.0 

14.2 

0.0 

136.9 

4320. 

3648. 

11270. 

9724. 

0.0 

82.4 

0.0 

65.9 

0.0 

10.3 

0.0 

91.1 

4880. 

4922. 

13050. 

13243. 

0.0 

132.8 

0.0 

101.9 

0.0 

11.8 

0.0 

136.9 

4710. 

4865. 

11850. 

9676. 

0.0 

86.4 

0.0 

73.4 

0.0 

8.6 

0.0 

89.8 

4000. 

5992. 

10180. 

12672. 

0.0 

130.5 

0.0 

112.5 

0.0 

10.7 

0.0 

134.1 

3760. 

4865. 

9340. 

9710. 

0.0 

87.2 

0.0 

74.3 

0.0 

8.7 

0.0 

89.8 

4960. 

5992. 

12120. 

13023. 

0.0 

130.2 

0.0 

111.4 

0.0 

10.6 

0.0 

134.1 

3590. 

2363. 

10790. 

7323. 

0.0 

89.5 

0.0 

77.1 

0.0 

18.6 

0.0 

91.6 

3670. 

2972. 

11760. 

10059. 

0.0 

134.3 

0.0 

112.6 

0.0 

21.6 

0.0 

138.3 
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Table 10. Continued 


gra ::< 

YIELD 

BIOMASS 

N UPTAKE 

GRAIN N 

NUPTAXE 

GRAIN 

PROTEIN 

AN - 4 '"S' 

3 N UPTAKE 

0 

P 

0 

P 

0 

P 

0 

P 

0 

P 

0 

P 

4070. 

2794. 

11050. 

9196. 

0.0 

84.0 

0.0 

67.2 

0.0 

13.7 

0.0 

91.6 

4440. 

3728. 

12170. 

12667. 

0.0 

134.1 

0.0 

110.5 

0.0 

16.9 

0.0 

138.3 

1220. 

4415. 

10700. 

10036. 

89.5 

106.3 

56.5 

92.9 

10.0 

12.0 

45.3 

109.7 

3940. 

4855. 

11701. 

12443. 

112.9 

149.5 

73.3 

113.3 

10.6 

13.3 

86.7 

134.6 

3940 . 

4963. 

7905. 

10480. 

63.4 

104.2 

58.8 

90.6 

8.5 

10.4 

75.3 

103.4 

4840. 

5954. 

11335. 

13761. 

104.9 

156.5 

88.3 

136.8 

10.4 

13.1 

93.1 

156.4 

3840. 

3348. 

9600. 

9451. 

83.7 

111.0 

67.4 

95.2 

10.0 

16.2 

70.7 

110.4 

4360. 

3837. 

10847. 

12030 . 

104.3 

152.9 

83.4 

131.9 

10.9 

19.6 

98.7 

157.4 

S040. 

4742. 

12642. 

12979. 

118.9 

144.1 

92.8 

96 : 

10.5 

11.6 

0.0 

0.0 

4030. 

4742. 

11060. 

12979 . 

112.6 

143.7 

78.5 

96.5 

11.1 

11.6 

0.0 

0.0 

3420. 

3934. 

9790. 

11247. 

98.3 

144.0 

68.4 

89.7 

11.4 

13.0 

0.0 

0.0 

1920. 

2524. 

4950. 

6366. 

41.9 

65.5 

31.0 

51.4 

9.2 

11.6 

0.0 

0.0 

4930. 

5620. 

12900. 

12059. 

127.1 

136.6 

97. 7 

119.3 

11.3 

12.1 

0.0 

0.0 

4080. 

5620. 

10124. 

12010. 

101.1 

133.9 

81.6 

117.3 

11.4 

11.9 

0.0 

0.0 

3470. 

5288. 

9582. 

11624 . 

90.6 

128.8 

68.2 

112.3 

11.2 

12.1 

0.0 

0.0 

1920. 

2987. 

4800. 

6310. 

39.7 

64.6 

31.7 

56.1 

9.4 

10.7 

0.0 

0.0 

4340 . 

3798. 

1C344 . 

12764 . 

95.1 

146.8 

78.1 

115.9 

9.8 

17.4 

0.0 

0.0 

4280. 

3798. 

9909. 

12764 . 

94.4 

146.4 

77.3 

115.9 

10.3 

17.4 

0.0 

0.0 

1150. 

3131. 

8600. 

10354. 

85.9 

143.9 

62.4 

98.3 

11.3 

17.9 

0.0 

0.0 

1340. 

1923. 

3968. 

5950. 

33.4 

65.9 

25.9 

55.3 

9.6 

16.4 

0.0 

0.0 

7520. 

5360. 

16460. 

11593 . 

179.6 

123.6 

143.8 

107.2 

10.9 

11.4 

166.5 

132.8 

8120. 

5828. 

18230. 

13883. 

198.6 

178.3 

156.7 

157.5 

11.0 

15.4 

212.6 

184.5 

7280. 

5876. 

17430. 

13940. 

237.6 

185.1 

178.8 

161.8 

14.0 

15.7 

0.0 

184.5 

7360. 

6652. 

16640. 

15856. 

218.0 

223.1 

169.2 

196.1 

13.1 

16.8 

225.6 

229.6 

6990. 

7018. 

16750. 

16481 . 

258.6 

2*4.3 

182.7 

206.8 

14.9 

16.8 

247.7 

252.2 

6570. 

5335. 

15870. 

11319. 

152.4 

125.1 

118.7 

110.4 

10.3 

11.8 

139.7 

125.4 

7750. 

6588. 

18180. 

14314. 

208.4 

177.1 

150.9 

157.2 

11.1 

13.6 

189.9 

180.6 

7650. 

6614. 

18180. 

14330. 

259.6 

181.4 

178.5 

161.3 

13.3 

13.9 

0.0 

180.6 

6910. 

7247. 

15380. 

15286. 

207.8 

201.9 

154.0 

181.8 

12.7 

14.3 

224.3 

204.9 

6340. 

7317. 

15380. 

15425. 

240.8 

206.0 

155.7 

184.9 

14.0 

14.4 

247.7 

209.1 

1100. 

692. 

2060. 

1366. 

27.5 

14.4 

25.5 

12.7 

13.2 

10.5 

0.0 

0.0 

2310. 

2881. 

4540. 

6049. 

48.3 

63.4 

40.9 

55.1 

10.1 

10.9 

0.0 

0.0 

3890- 

4221. 

7170- 

9365. 

90.0 

103.4 

77.8 

89.6 

11.4 

12.1 

0.0 

0.0 

4850. 

4332 . 

8830. 

9790. 

118.0 

114.3 

105.5 

100.3 

12.4 

13.2 

0.0 

0.0 

3530. 

3781. 

6730. 

8182. 

90.3 

89.3 

76.2 

77.6 

12.3 

11.7 

0.0 

0.0 

3550. 

3781. 

6690. 

8155. 

95.2 

87.9 

74.7 

76.3 

12.0 

11.5 

0.0 

0.0 

2960. 

2990. 

5800. 

6137. 

87.5 

67.3 

69.1 

58.8 

13.3 

11.2 

0.0 

0.0 

2990. 

2990. 

5930. 

6137. 

92.5 

67.4 

77.6 

59.3 

14.8 

11.3 

0.0 

0.0 

4160. 

4221. 

7860. 

9368. 

105 9 

103.8 

84.7 

90.3 

11.6 

12.2 

0.0 

0.0 

6290. 

8171. 

15720. 

16500 . 

147.0 

186.4 

108.1 

166.3 

9.8 

11.6 

90.0 

187.0 

6910. 

8228. 

15710. 

16578. 

183.0 

187.9 

138.2 

167.4 

11.4 

11.6 

100.0 

187.7 

8560. 

8229. 

18210. 

16579. 

226.0 

187.9 

171.2 

167.5 

11.4 

11.6 

130.0 

188.9 

8580. 

8229. 

17870. 

16579. 

248.0 

187.9 

192.7 

167.5 

12.8 

11.6 

142.0 

188.9 

4170. 

6132. 

9100. 

11896. 

71.0 

114.7 

54.1 

100.0 

7.4 

9.3 

0.0 

0.0 

5640. 

7043. 

12200. 

13781. 

105.0 

140.9 

81.1 

123.6 

8.2 

10.0 

0.0 

0.0 

6020. 

7933. 

13000. 

15857. 

126.0 

175.5 

94.0 

153.1 

8.9 

11.0 

0.0 

0.0 

6260. 

8297. 

13000. 

17169. 

162.0 

207.3 

114.2 

165.9 

10.4 

11.4 

0.0 

0.0 

4620. 

6949. 

10500. 

13363. 

89.0 

135.1 

66.5 

118.3 

8.2 

9.7 

0.0 

0.0 


43 



T. bl« 10. Continued. 


GRAIN 

YIELD 

BIOMASS 

N UPTAKE 

GRAIN N 

NUPTAKE 

CRAIN 

PROTEIN 

ANTHESIS N UPTAKE 

0 

P 

0 

P 

0 

P 

0 

P 

0 

P 

0 

P 

$870. 

7626. 

13000. 

15188. 

118.0 

162.9 

87.5 

141.8 

8.5 

10.6 

0.0 

0.0 

$230. 

8167. 

13400. 

16715. 

147.0 

195.1 

104.9 

161.9 

9.6 

11.3 

0.0 

0.0 

6360. 

8511. 

1’ 00. 

17773. 

166.0 

224.9 

121.6 

170.2 

10.9 

11.4 

0.0 

0.0 

1400. 

2975. 

8680. 

9343. 

74.3 

124.1 

56.7 

84.0 

9.5 

16.1 

81.9 

129. 3 

3730. 

3022 

9190. 

9854. 

76.0 

134.6 

69.4 

88.5 

10.6 

16.7 

79.5 

140.0 

3460. 

3C«.3. 

8910. 

10206. 

86.3 

135.8 

66.2 

91.8 

10.9 

17.3 

88.7 

147.9 

3410. 

3' 17. 

7920. 

10353 . 

68.4 

123.4 

60.4 

106.5 

10.1 

15.5 

70.6 

128.3 

1660 

3: 97. 

8690. 

10857 . 

74.1 

134.3 

68.7 

116.2 

10.7 

17.0 

77.8 

138.3 

3650. 

41*0. 

8550. 

11287 . 

92.2 

141.8 

69.8 

122.6 

10.9 

17.3 

70.5 

146.9 

5810. 

42b . 

14470. 

13888. 

121.7 

211.0 

97.9 

161.5 

9.6 

21.5 

129.7 

220.6 

5600. 

4377. 

14370. 

14410. 

133.2 

222.8 

101.2 

165.1 

10.3 

21.5 

141.8 

237.9 

5490. 

4301. 

14670. 

14876. 

158.7 

238.8 

121.4 

169.0 

12.6 

22.4 

185.1 

254.2 

5530. 

4414. 

14770. 

15406. 

179.6 

257.8 

134.9 

173.5 

13.9 

22.4 

197.9 

272.2 

5190. 

6532. 

12660. 

16305. 

111.5 

213.3 

89.2 

186.8 

9.8 

16.3 

124.8 

221.4 

5690. 

6976. 

14100. 

17216. 

137.9 

232.0 

107.8 

203.2 

10.8 

16.6 

139.5 

240.7 

5310. 

6897. 

13780. 

17835. 

159.8 

251.6 

118.3 

216.6 

12.7 

17.9 

145.3 

260.7 

5190. 

7117. 

13600. 

18460. 

158.1 

271.2 

122.0 

224.7 

13.4 

18.0 

161.0 

280.8 

2070. 

2252. 

5694. 

3973. 

71.0 

46.6 

50.8 

42.3 

14.0 

10.7 

0.0 

0.0 

3772. 

4409. 

8412. 

8349. 

129.7 

95.3 

100.6 

85.1 

15.2 

11.0 

0.0 

0.0 

3450. 

5057. 

8459. 

9469. 

116.7 

107.6 

87.2 

95.8 

14.4 

10.8 

0.0 

0.0 

1470. 

4899. 

8736. 

9074. 

127.4 

104.7 

93.8 

93.7 

15.4 

10.9 

0.0 

0.0 

3631. 

5069. 

9292. 

9481. 

136.9 

109.1 

100.0 

96.9 

15.7 

10.9 

0.0 

0.0 

804. 

313. 

2821. 

543. 

32.8 

6.6 

20.7 

6.0 

14.7 

11.0 

0.0 

0.0 

992. 

1141. 

4286. 

1959. 

50.5 

24.6 

25.8 

22.4 

14.8 

11-2 

0.0 

0.0 

1159. 

1558. 

8241. 

2682. 

77.5 

33.8 

31.5 

30.6 

15.5 

11.2 

0.0 

0.0 

1240. 

1727. 

8322. 

2976. 

90.8 

37.7 

34.2 

34.2 

15.7 

11.3 

0.0 

0.0 

1521. 

777. 

3829. 

1449. 

53.7 

17.5 

39.0 

15.8 

14.6 

11.6 

0.0 

0.0 

1307. 

1067. 

4086. 

2001. 

65.1 

24.0 

37.6 

21.5 

16.4 

11-5 

0.0 

o.c 

978. 

1161. 

3264. 

2178. 

58.3 

26.1 

28.3 

23.4 

16.5 

11.5 

0.0 

0.0 

1005. 

1258. 

34C3. 

2270. 

60.9 

27.1 

29.1 

24.5 

16.5 

11.1 

0.0 

0.0 

2352. 

3437. 

0 . 

9420. 

48.2 

88.9 

35.5 

62.7 

8.6 

10.4 

0.0 

90.6 

3024. 

4090. 

0 . 

12657. 

73.9 

139.8 

48.3 

74.6 

9.1 

10.4 

0.0 

143-5 

3696. 

4202. 

0 . 

13874. 

96.3 

169.4 

67.4 

76.7 

10.4 

10.4 

0.0 

174.4 

3584. 

4375. 

0 . 

14936. 

132.2 

204.3 

91.2 

79.8 

14.5 

10.4 

0.0 

209.7 

2240. 

2978. 

0 . 

10871. 

50.4 

131.0 

33.8 

56.4 

8.6 

10.8 

0.0 

133.1 

2800. 

3145. 

0 . 

12223. 

65.0 

163.4 

44.7 

59.6 

9.1 

10.8 

0.0 

16* . 3 

2688. 

3254. 

0 . 

13201. 

67.2 

184.8 

49.0 

61.7 

10.4 

10.8 

0.0 

195.7 

2352. 

3375. 

0 . 

14663. 

118.7 

235.9 

59.8 

63.9 

14.5 

10.8 

0.0 

243.2 

670. 

447. 

2016. 

869. 

14.6 

9.6 

0.0 

8.5 

0.0 

10.9 

0.0 

0.0 

1337. 

2041. 

4256. 

4140. 

31.3 

47.8 

0.0 

42.6 

0.0 

11.9 

0.0 

0.0 

2866. 

2101. 

10416. 

10035 . 

154.5 

178.7 

95.5 

53.4 

19.0 

14.5 

222.9 

186.3 

2240. 

2209. 

9744. 

10996 . 

155.7 

214.1 

60.5 

58.1 

15.4 

15.0 

227.4 

222.1 

2128. 

2281. 

9632. 

11332. 

150.1 

234.4 

57.9 

60.0 

15.5 

15.0 

198.2 

242.6 

2016. 

2273. 

9408. 

11313. 

151.2 

241.2 

70.7 

59.8 

20.0 

15.0 

212.8 

249.2 

2800. 

2619. 

10864. 

11367. 

154.5 

178.6 

94.3 

66.6 

19.2 

14.5 

222.9 

190.2 

2576. 

2820. 

10528. 

12342. 

155.7 

214.3 

85.0 

71.7 

18.8 

14.5 

227.4 

223.2 

2352. 

2891. 

10080. 

12670. 

150.1 

234.6 

78.4 

73.5 

19.0 

14.5 

198.2 

243.4 

2240. 

2878. 

10080. 

12620 . 

151.1 

241.3 

78.6 

73.2 

20.0 

14.5 

212.8 

251. 1 
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Table 10. Continued 


GRAIN 

YIELD 

BIOMASS 

N UPTAKE 

GRAIN N 

NUPTAKE 

GRAIN 

PROTEIN 

ANTHESIS N UPTAKE 

0 

P 

0 

P 

0 

P 

0 

P 

0 

P 

0 

P 

2866. 

2108. 

10416. 

10087 . 

154.5 

184.6 

95.5 

53.6 

19.0 

14.5 

222.9 

192.3 

2260. 

2204. 

9744. 

10997 . 

155.7 

218.3 

60.5 

58.0 

15.4 

15.0 

227.4 

226.4 

2128. 

2275. 

9632. 

11318. 

150.1 

237.6 

57.9 

59.9 

15.5 

15.0 

198.2 

245.6 

2016. 

2265. 

9408. 

11297. 

151.2 

243.1 

70.7 

59.6 

20.0 

15.0 

212.8 

251.5 

2800. 

2621. 

10864. 

11406. 

154.5 

184.8 

94.3 

66.7 

19.2 

14.5 

222.9 

195.6 

2576. 

2815. 

10528. 

12327 . 

155.7 

218.5 

85.0 

71.6 

18.8 

14.5 

227.4 

227.4 

2352. 

2884. 

10080. 

12647 . 

150.1 

237.7 

78.4 

73.4 

19.0 

14.5 

198.2 

246.9 

2240. 

2869. 

10080. 

12605. 

151.1 

243.3 

78.6 

73.0 

20.0 

14.5 

212.8 

252.6 

329. 

139. 

0 . 

265. 

8.3 

2.3 

0.0 

1.9 

0.0 

7.8 

0.0 

0.0 

531. 

1043. 

0 . 

2077. 

16.2 

15.6 

0.0 

12.8 

0.0 

7.0 

0.0 

0.0 

766. 

888. 

0 . 

1749. 

25.2 

13.2 

0.0 

10.9 

0.0 

7.0 

0.0 

0.0 

1082. 

1270. 

0 . 

2492. 

27.4 

19.2 

0.0 

16.0 

0.0 

7.2 

0.0 

0.0 

1344. 

1425. 

0 . 

2642. 

30.8 

20.8 

0.0 

17.8 

0.0 

7.1 

0.0 

0.0 

3400. 

3787. 

5650. 

7488. 

49.0 

60.0 

47.1 

50.5 

7.9 

7.6 

40.0 

65.9 

4000. 

4943. 

7550. 

9904. 

57.0 

87.3 

51.9 

74.6 

7.4 

8.6 

40.0 

90.8 

5400. 

5803. 

10000. 

11845. 

83.0 

108.2 

76.7 

92.6 

8.1 

9.1 

0.0 

115.6 

6400. 

6398. 

12430. 

13576. 

107.0 

133.7 

96.6 

115.6 

8.6 

10.3 

96.0 

137.6 

6600. 

7162. 

13250. 

15281 . 

123.0 

160.0 

110.0 

138.2 

9.5 

11.0 

96.0 

163.8 

7000. 

7754. 

14250. 

16608. 

153.0 

179.6 

136.3 

156.4 

11.1 

11.5 

0.0 

185.0 

6600. 

8018. 

13000. 

17629. 

149.0 

207.6 

134.3 

177.2 

11.6 

12.6 

172.0 

211.8 

7000. 

8105. 

14060. 

17751. 

197.0 

220.7 

158.4 

182.0 

12.9 

12.8 

172.0 

225.2 

3600. 

3788. 

7300. 

8039. 

62.0 

60.6 

59.4 

47.2 

9.4 

7.1 

40.0 

67.0 

4200. 

4739. 

8950. 

10355. 

69.0 

89.7 

64.1 

67.3 

8.7 

8.1 

40.0 

92.2 

4800. 

5409. 

10000. 

12066. 

72.0 

113.1 

67.4 

76.9 

8.0 

8.1 

0.0 

117.4 

5800. 

5821. 

12700. 

13314. 

98.0 

134.3 

87.5 

82.7 

8.6 

8.1 

96.0 

138.2 

6100. 

6209. 

13000. 

14457. 

111.0 

156.3 

98.5 

88.2 

9.2 

8.1 

96.0 

161.6 

6800. 

6610. 

15050. 

15564. 

153.0 

181.9 

126.5 

93.9 

10.6 

8.1 

0.0 

186.6 

6200. 

6991. 

1 4300. 

16582. 

150.0 

207.4 

121.8 

99.3 

11.2 

8.1 

172.0 

212.1 

6500. 

7111. 

14500. 

17384. 

172.0 

230.4 

138.0 

103.5 

12.1 

8.3 

172.0 
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Summary 


CERES -Wheat -N is designed to be used as a management - 
oriented model for a diversity of applications in many 
environments. To make the model useful for such a wide 
audience, the inputs must be minimal and they must be 
reasonably easy to attain or estimate from standard agricul- 
tural experimental practice. Given these constraints, the 
model simulates crop growth and response to fertilizer 
reasonably reliably. The rigorous statistical analysis 
indicated that some problems have yet to be resolved with the 
prediction of N uptake and grain protein concentration in 
some data sets. Some tendency to overpredict biomass and N 
uptake early in the season was also noted. Further testing 
and refinement of this area of the model may be beneficial. 
Further data sets are required to test more reliably the soil 
N components of the model. 
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Figure B.5. Comparison of predicted and observed grain protein response 
to applied N in individual data sets. 
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Figure B.6. Compariscn of predicted and observed N uptake response to 
applied N in individual data sets. 
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Figure B.8. Comparison of predicted 
uptake for individual data sets. 
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